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PREFACE 


This  report  was  prepared  by  the  Lockheed-California  Company,  Burbank, 
California  for  the  Structural  Integrity  Branch,  Structures  Division,  Air 
Force  Flight  Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio 
under  Contract  F33615-75-C-3087.  The  work  described  herein  is  a continu- 
ing part  of  the  Air  Force  Systems  Command's  exploratory  development  program 
to  investigate  all  practical  load  combinations  affecting  crack  growth,  to 
meet  the  damage  tolerance  requirements  for  safety  of  flight  aircraft  struc- 
tures. The  work  was  directed  under  Project  FI  4 56- 7 5-574/1*8611  and  task  486U02. 
Mr.  D.  Simpkins  of  the  Structural  Integrity  Branch  was  the  task  engineer. 

This  report  concludes  the  work  on  Contract  F33615-75-C-3087,  which  covered 
a period  from  June,  1975  to  July  1976. 

The  development  of  the  flight-by-flight  load  spectra  and  the  analysis  of 
the  flight  load  crack  growth  was  performed  by  Dr.  S.T.  Chiu,  the  acoustic 
testing  and  data  reduction  by  Mr.  J.  L.  Hayward  and  the  crack  propagation 
testing  and  data  reduction  by  Mr.  C.J.  Looper.  Mr.  J.  Soovere  was  the 
project  leader. 


The  manuscript  was  released  by  the  authors  for  publication  on  29th  July  1976. 
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over-bar  denotes  the  root  mean  squa.re  value  or  constants 
associated  with  the  root  mean  square  stress  except  in  the 
above  list  of  symbols 

evaluated  at  given  constants  of  x and  y or  any  other 
variable  so  defined 


Subscript 


A 

denotes 

acoustic 

eff 

denotes 

effective 

j 

denotes 

cycle  number 

M 

denotes 

mass  attached 

max 

denote  s 

maximum  value 

of 

min 

denotes 

minimum  value 

of 

OL 

denotes 

overload 

S 

denotes 

shaker 
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I. 


INTRODUCTION 


Crack  propagation  studies  are  employed  in  the  aircraft  industry 
to  formulate  inspection  and  maintenance  schedules  for  "in-service"  monitor- 
ing to  meet  the  damage  tolerance  requirements  (Reference  l)  for  safety  of 
flight  structures.  Interest  is,  therefore,  centered  on  all  practical  load- 
ing combinations,  such  as  the  combination  of  flight-by- flight  loads  with 
acoustic  loads,  which  are  not  accounted  for  in  current  crack  propagation 
analysis  and  which  may  have  a significant  effect  on  the  crack  growth  behavior. 
The  potential  effect  of  the  combined  loading  has  assumed  a greater  signi- 
ficance with  the  advent  of  the  Advanced  Medium  Transport  Aircraft.  In  these 
aircraft,  due  to  the  close  proximity  of  the  jet  engine  exhaust  to  the  air- 
craft structure,  a significant  level  of  acoustic  excitation  of  a long  dura- 
tion will  be  encountered  during  flight  in  addition  to  the  high  level  take- 
off and  landing  noise  levels  which  normally  impose  the  sonic  fatigue  design 
requirement  for  critical  structures. 

Crack  propagation  in  aircraft  structures  due  to  either  flight  loads 
or  acoustic  loads  has  been  studied  extensively,  but  independently,  over  the 
past  two  decades.  The  separated  development  is  due  in  part  to  the  analysis 
methods  which  lean  heavily  on  their  respective  fatigue  analysis  procedures. 

The  flight  loads  can  introduce  significant  stress  levels  in  both  heavy  and 
light  load  carrying  structures  whereas  significant  stress  levels  can  only 
be  produced  in  relatively  light  panel  type  structures  by  acoustic  loading. 

The  loading  action  is  also  different  for  the  acoustic  and  the  flight  loads. 
The  flight  loads  introduce  predominantly  inplane  tension  and  compression 
stresses  whereas  the  acoustic  loading  introduces  predominantly  bending 
stresses,  in  the  local  structure. 

In  the  study  of  crack  propagation  due  to  flight  loads,  the  re- 
presentation of  the  flight  loads  has  progressed  from  constant  amplitude 
loading  (Reference  2)  to  the  more  complex  flight-by- flight  loading  (Refer- 
ences 3 to  6).  To  improve  the  accuracy  of  the  latter  analyses,  methods 
(References  7,  8 and  9)  which  include  the  effects  of  crack  retardation  due 
to  infrequent  high  loads  have  been  developed.  The  test  specimen  configura- 
tions have  progressed  from  the  simple  center  crack  specimen  to  the  more 
complex  structural  joint  specimen  with  precracks  in  the  fastener  holes 
(References  10  and  ll). 

The  study  of  crack  propagation  due  to  acoustic  loading  was  initiat- 
ed (Reference  12)  in  connection  with  turbulent  boundary  layer  excitation  of 
the  fuselage  panels  in  a supersonic  transport  aircraft  at  high  dynamic 
pressures.  The  turbulent  boundary  layer  was  simulated  by  acoustic  loading, 
with  the  constant  fuselage  pressurization  stresses  represented  by  simul- 
taneously applied  uniaxial  (References  12  and  13)  and  later  (Reference  l4) 
by  biaxially  inplane  loads.  Again,  test  panels  with  a center  crack  normal 
to  the  principal  acoustic  and  inplane  hoop  loading  were  employed. 
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Subsequent  studies  employing  the  more  realistic  panel  edge  cracks,  indicat- 
ed that  significant  crack  growth  can  be  obtained  due  to  peak  acoustic  load- 
ing during  takeoff  in  sequence  with  the  ground- air- ground  cycle  (Reference 
15)  and  indeed  from  the  takeoff  acoustics  alone  (Reference  1 6).  The  basic 
analysis  procedure  (Reference  l6)  is  based  on  random  response  theory  employ- 
ing the  Rayleigh  distribution  of  stress  peaks  associated  with  a predominant 
single  mode  response.  The  test  data  indicates  that  loading  frequency  and 
crack  closure  affect  the  crack  growth  rate. 

When  flight  loads  and  acoustic  loads  are  applied  simultaneously, 
interaction  effects  can  occur  which  may  affect  the  crack  growth  rate.  This 
program  represents  a preliminary  investigation  into  the  study  on  the  effect 
of  combined  acoustic  and  flight  loads  on  crack  growth. 

The  objectives  of  this  program  are  to  determine 

o if  combined  acoustic  and  flight  loads  have  a signi- 
ficant effect  on  crack  growth 

o whether  any  significant  effects  can  be  predicted  by 
current  crack  propagation  analysis  methods 

As  no  test  facilities  are  available  in  which  flight  and  acoustic  loads 
can  be  simultaneously  applied,  electromagnetic  shaker  simulation  of  acoustic 
loading  is  employed  in  this  program.  Other  requirements  imposed  on  this  pro- 
gram are  the  need  for  real  time  simulation  of  the  duration  of  the  acoustic 
excitation  as  far  as  practical  and  the  application  of  significant  compression 
loads  and  tension  flight  loads  on  the  test  specimen.  Airplane  skin  bay  type 
structures,  which  are  subjected  to  significant  compression  loads,  are  gener- 
ally loaded  axially  with  the  smallest  panel  dimension  normal  to  the  applied 
load  to  obtain  maximum  resistance  to  panel  buckling.  With  the  simultaneous 
application  of  acoustic  and  axial  flight  loading,  a predominantly  biaxial 
stress  field  is  obtained  at  the  center  of  the  skin  bay  and  at  the  center  of 
the  longest  panel  side.  At  these  locations,  the  axial  flight  load  stresses 
are  normal  to  the  predominant  acoustic  load  stresses.  A greater  variability 
and  higher  stress  levels  are  usually  encountered  with  the  flight -by-flight 
loading  in  comparison  to  the  acoustic  loading.  Consequently,  the  programmed 
crack  orientation  used  in  the  combined  loads  tests  in  this  program,  is  selected 
on  the  premise  that  the  flight  loading  is  the  more  critical  loading  condition 
and,  therefore,  more  likely  to  initiate  the  crack. 

This  report  on  the  research  program  contains  the  specimen  design 
considerations  including  panel  response  theory  and  the  impact  of  electro- 
magnetic shaker  simulation  of  the  acoustic  loading  on  the  panel  design  in 
Section  II;  the  development  of  the  loading  time  history  based  on  truncated 
runway  and  gust  spectra  in  Section  III;  a summary  of  current  crack  propagation 
prediction  methods  and  crack  growth  calculations  in  Section  IV;  the  test  pro- 
gram in  Section  V;  the  correlation  and  evaluation  of  the  test  and  analytical 
results  in  Section  VI;  and,  finally,  the  conclusions  in  Section  VII. 
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II  DESIGN  CONSIDERATIONS 


1.  Introduction 


The  purpose  of  the  program  is  to  study  crack  growth  behaviors  in  simu- 
lated aircraft  structures  due  to  combined  acoustic  and  flight  loads.  The 
design  of  the  test  specimen  must,  therefore,  be  such  that  significant 
stresses  can  be  achieved  in  the  region  of  the  programmed  crack  by  both  load- 
ing actions.  The  test  specimen  must,  in  addition,  withstand  compression 
loads  without  buckling. 

These  requirements  can  be  satisfied  by  a single  panel  test  specimen  sup- 
ported by  stiffeners  on  both  sides.  Since  an  electromagnetic  shaker  is  used 
to  simulate  the  acoustic  excitation,  the  shaker  coil  and  the  drive  rod 
weight  will  significantly  affect  the  panel  response  frequency  and  stresses 
and,  consequently,  the  panel  dimensions.  The  panel  dimensions  are  equally 
affected  by  the  compression  loads.  A theoretical  study  of  the  panel  response 
is,  therefore,  necessary  to  define  the  dimensions  of  the  test  panel. 

2.  Panel  Frequency 

The  dynamic  response  of  stiffened  panel  type  aircraft  structures  represents 
a complex  analytical  problem  due  to  the  flexibility  of  the  panel  edges.  The 
fundamental  frequency  of  the  stiffened  panels  falls  halfway  between 
the  fixed-fixed  and  the  simply  supported  panel  frequencies,  tending  towards 
the  corresponding  simply  supported  frequency  for  progressively  higher  panel 
modes.  When  typical  stiffened  aircraft  panels  are  excited  by  jet  noise,  pre- 
dominant single  mode  response  is  observed  (References  17  and  l8).  In  acous- 
tic progressive  wave  tunnels,  used  to  simulate  jet  noise  excitation,  the 
second  axial  panel  mode  can  also  be  excited  if  the  panel  longest  side  is 
parallel  to  the  tunnel  axis.  In  general,  the  use  of  only  the  fundamental  mode 
in  panel  response  studies  of  random  acoustic  excitation  is  adequate. 

The  expression  for  the  simply  supported  panel  fundamental  frequency,  when 
modified  by  an  empirical  factor  to  account  for  the  edge  flexibility 
(Reference  19),  can  be  used  to  predict  the  stiffened  panel  fundamental  fre- 
quency. The  advantages  are  its  simplicity  and  accuracy.  This  method  will 
now  be  extended  to  include  the  effects  of  an  axial  inplane  load  and  the 
addition  of  a panel  center  mass,  representing  the  weight  of  the  shaker  coil 
and  drive  rod,  on  the  panel  fundamental  frequency.  For  maximum  compression 
load  carrying  capability,  the  smallest  panel  dimension  must  be  normal  to 
the  direction  of  the  applied  compression  load. 

2.1  Effect  of  Inplane  Load 

The  differential  equation  for  the  deflection,  W,  of  a flat  plate  sub- 
ject to  axial  inplane  load,  Nx,  per  unit  width  is  given  by  (References 
20  and  21,  pp.  276) 
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where 
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D = - E h3 
12(1-  v 2 ) 


is  the  flexural  rigidity  of  the  panel,  q is  the  external  transverse  surface 
force  per  unit  area,  h is  the  panel  thickness,  E is  the  Young's  modulus 
and  v is  the  Poisson's  ratio.  For  dynamic  response,  q is  replaced  by  the 
inertia  force  given  by 


where  p is  the  panel  density.  The  panel  coordinates  are  defined  in 
Figure  1.  The  deflection  in  any  one  mode  of  a simply  supported  panel  is 


is  the  simply  supported  panel  mode  shape,  m and  n are  the  axial  and  lateral 
mode  number,  a is  the  panel  length,  and  b is  the  panel  width.  For  simple 
harmonic  motion,  £(t),is  given  in  terms  of  the  generalized  coordinate, f , 

■by 


Substitute  equations  2 to  5 into  equation  1.  After  the  necessary 
algebraic  manipulation,  the  following  expression  is  obtained  for  the 
panel  frequency 


With  zero  inplane  load,  equation  6 reduces  to  the  expression  for  the 
frequency  of  a simply  supported  panel  given  by 
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W = w(x,y ) *£(t) 
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where 
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5(t)  - £ sin  wt 
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FIGURE  1.  PANEL  COORDINATES 


Substituting  equation  7 back  into  equation  6 results  in 
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The  critical  panel  compressive  stress,  axc  > is  obtained  when  the  panel 
frequency  becomes  zero.  Then 


xc 
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which  is  the  relationship  between  the  unloaded  panel  frequency  and  the  criti- 
cal axial  buckling  stress. 


2.2  Effect  of  Panel  Edge  Flexibility 

The  effect  of  edge  flexibility  on  panel  frequency  is  accounted  for  by  the 
introduction  of  an  empirical  constant  K (Reference  19)  into  equation  7 as 
illustrated  below 


fr 
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For  the  fundamental  panel  mode  K is  approximately  equal  to  1.48.  On  re- 
placing the  simply  supported  panel  frequency  in  equation  8 with  the  stiffen- 
ed panel  frequency,  equation  10,  and  on  rearranging  the  terms,  the  expres- 
sion for  the  stiffened  panel  frequency  with  inplane  load  is  obtained  and 


given  by 


f = f„ 
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The  critical  buckling  stress  is  a function  of  panel  aspect  ratio.  The 
lowest  critical  buckling  stress  may  occur  with  an  axial  panel  mode  higher 
than  the  fundamental.  Only  limited  data  (Reference  19)  is  available  on  the 
higher  order  modal  frequencies  for  typical  aircraft  panels.  This  data  is 
insufficient  for  constructing  a frequency  correction  curve  for  the  higher 
order  axial  panel  modes.  The  frequency  correction  factor,  for  the  higher 
order  axial  panel  modes, can,  however,  be  obtained  from  the  theory  in  Refer- 
ence 21,  page  117,  which  predicts  the  effect  of  panel  edge  rotational  re- 
straint on  panel  frequency.  The  two  longer  panel  edges  are  rotationally 
restrained  whereas  the  two  shorter  edges  are  simply  supported.  A represen- 
tative rotational  edge  stiffness  parameter  can  be  •extracted  from  the  theory 
by  requiring  the  ratio  of  the  stiffened  panel  frequency  to  the  simply  sup- 
ported panel  frequency  to  equal  1.48  at  a panel  aspect  ratio  of  two.  There- 
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after,  the  remainder  of  the  frequency  correction  curve  is  readily  obtained 
from  the  theory  and  is  illustrated  in  Figure  2 as  a function  of  effective 
aspect  ratio  defined  by  (a/mb).  Limited  test  data  is  included  in  Figure  2 
for  comparison. 


2.3  Effect  of  Added  Panel  Center  Mass 


The  effect  of  an  added  panel  center  mass  M on  the  panel  frequency  can  be 
computed  considering  the  generalized  mass  of  the  panel  and  the  added  weight. 
The  generalized  mass  Mr  for  the  panel  is  given  by 


M 
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The  generalized  mass  of  the  added  mass  is  equal  to  M for  a simply  supported 
mode  shape  as  the  modal  displacement  is  unity  at  the  center  of  the  panel. 
The  corresponding  panel  generalized  mass  is  equal  to  0.25  pabh.  The  ratio 
of  the  generalized  mass  of  the  panel  to  the  generalized  mass  of  the  added 
mass  and  the  panel  is 


pabh 

pabh  + 4M 


The  relationship  between  the  simply  supported  panel  frequency  with  (Refer- 
ence 21,  pp.  1V7)  and  without  the  attached  mass  follows  directly  and  is 
given  by 


M 


= f 


/ pabh  \ 
s vpabh  + hul 
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2.b  Panel  Frequency  Design  Equation 


The  addition  of  the  center  mass  will  not  affect  the  critical  buckling  stress 
which  is  still  related  to  the  stiffened  panel  frequency  without  the  attached 
mass.  The  final  expression  for  the  stiffened  panel  frequency  with  both  the 
axial  inplane  load  and  the  attached  mass  is  given  by 


where  fr  is  given  by  equation  10. 
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3.  Panel  Stress 

In  typical  stiffened  aircraft  panels  which  are  subjected  to  acoustic  load- 
ing, the  greatest  panel  bending  stresses  occur  in  the  fundamental  mode  in  a 


7 


On 


8 


FIGURE  2.  FREQUENCY  CORRECTION  FACTOR  FOR  HIGHER  ORDER  AXIAL  PANEL  MODES 


direction  normal  to  the  panel  edge  both  at  the  center  of  the  longest  panel 
edge  and  at  the  center  of  the  panel.  The  resulting  lateral  panel  stress 
distribution  is  generally  symmetrical  about  the  panel  centerline.  After 
peaking  at  the  panel  edge,  the  stress  distribution  falls  off  to  zero  some- 
where between  the  panel  edge  and  center  and,  thereafter,  to  a negative 
peak  at  the  panel  center.  The  magnitude  of  the  panel  edge  and  center 
stresses  are  comparable  for  the  level  of  edge  flexibility  encountered  in 
stiffened  aircraft  panel  structures.  The  panel  edge  is,  however,  more 
critical  for  fatigue  crack  initiation  because  of  the  stress  concentration 
effect  of  the  frame  edge  and  the  fastener  holes.  Sonic  fatigue  and  flight 
load  cracks  usually  emanate  from  fastener  holes  in  this  region  which  has 
been  chosen  as  the  location  for  the  programmed  cracks  in  this  study.  Con- 
sequently, attention  will  be  focused  on  panel  stresses  in  the  panel  edge 
region . 

Accurate  prediction  of  panel  edge  stresses  in  these  stiffened  panels  due  to 
acoustic  loading  is  difficult  as  the  magnitude  and  distribution  of  stresses 
are  affected  by  such  parameters  as  the  stiffener  geometry,  the  stiffener 
width  and  the  method  of  attachment  (References  23  and  2b).  Sufficient  empi- 
rical data  is  available  on  the  panel  edge  stresses  for  typical  stiffened 
panels  when  excited  by  random  acoustic  loading  (Reference  19) • Consequently, 
it  is  only  necessary  to  predict  changes  in  the  panel  edge  stress  due  to 
the  attached  mass,  damping  and  inplane  load  in  order  to  use  the  above  data 
for  design  purposes.  This  method  is  based  on  the  assumption  that  the  addi- 
tion of  the  center  mass  and  the  application  of  inplane  load  do  not  affect 
the  mode  shape  in  the  vicinity  of  the  panel  edge. 

The  addition  of  a panel  center  mass  and  the  use  of  the  shaker  to  simulate 
the  acoustic  excitation,  will  increase  the  force  requirement,  not  only  be- 
cause of  the  increased  inertia  force,  but  also  due  to  the  added  back  electro- 
motive force  (EMF)  induced  damping  from  the  shaker.  The  back  EMF  damping 
can  be  eliminated  by  switching  off  the  shaker  and  maintaining  open  Circuits 
during  the  acoustic  excitation  calibration  tests.  Thus,  the  combined  as 
well  as  independent  effects  of  added  mass  and  damping  must  be  determined. 

When  the  psuedo  random  flight  loads  are  simultaneously  applied  with  the 
simulated  random  acoustic  loading  a biaxial  stress  field  will  be  obtained 
in  the  panel.  Insight  into  the  nature  of  the  biaxial  stress  field  at  the 
panel  edge  can  be  obtained  by  considering  the  instantaneous  stress  and 
strain  state  within  the  panel. 

3.1  Panel  Stress  from  Acoustic  Excitation 

The  panel  edge  axial  and  lateral  modal  strains,  txe  and  «ye,  respectively, 
are  related  to  the  panel  mode  shape  w(x,y)  by 
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The  random  acoustic  spectrum  represented  by  G ( £1  , I'?  , w ) in  its  most 
general  form,  includes  information  on  the  distribution  of  the  acoustic 
pressure  field  over  both  the  surface  of  the  panel  and  the  frequency  range. 

| ^ and  | „ are  separation  distances  on  the  surface  of  the  panel  given  by 
(x2  _x]_)  (yg  -y^) , respectively,  and  w is  the  circular  frequency.  Using  the 
normal  modes  approach,  the  mean  square  panel  edge  strain  due  to  random 

acoustic  excitation  (References  25  and  2 6)  is  given  by 
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and  ur  is  the  panel  resonant  circular  frequency  and  6 is  the  viscous 
damping  coefficient.  The  double  area  integral  in  equation  17  represents 
the  generalized  load  or  force  spectrum.  As  the  white  noise  spectrum  level 
is  independent  of  frequency,  equation  17  can  readily  be  integrated  with 
respect  to  the  circular  frequency  (Reference  26),  using  contour  integration, 
to  obtain  the  mean  square  panel  edge  strain.  After  integration,  the  mean 
square  strain  is  given  by 


where 
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is  the  generalized  white  noise  load  spectrum 
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The  panel  edge  frames  are  assumed  to  be  stiff  in  bending.  Consequently, 
the  axial  strain  at  the  panel  edge  is  zero.  The  modal  panel  edge  stress 
and  strain  at  the  panel  edge  are  related  by 
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The  mean  square  panel  edge  stress  is,  therefore,  given  by 
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3.2  Effect  of  Added  Mass  on  Panel  Stress 


The  effect  of  the  added  panel  center  mass  is  to  reduce  the  panel  frequency 
•while  increasing  its  generalized  mass.  The  resulting  mean  square  panel  edge 
stress  normal  to  the  panel  edge  becomes 
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Thus,  the  ratio  of  mean  square  panel  edge  stress  with  the  added  mass  and 
without  the  added  mass  is 
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By  substituting  the  panel  frequency,  f<j,  at  zero  inplane  load  (equation  l4) 
into  the  above  equation,  it  becomes 


11 


- 2 

°yeM 

- 2 

a 

ye 


27 


Depending  on  the  magnitude  of  the  added  mass,  the  difference  in  the  panel 
frequencies  with  and  without  the  attached  mass  could  be  very  large.  A large 
difference  in  the  panel  frequency  could  affect  the  magnitude  of  both  the 
generalized  load  spectrum  and  the  damping  coefficient  (Reference  27)  as  both 
are  a function  of  frequency.  The  empirical  data  in  Reference  19  is  in  the 
form  of  design  charts  which  are  derived  on  the  assumption  of  a fully  correla- 
ted pressure  field.  Thus,  when  using  these  empirical  data  in  equation  27, 
the  generalized  load  spectra  GF^  and  GF  can  be  considered  to  be  equal  for 
the  same  acoustic  spectrum  level.  The  magnitude  of  the  panel  damping  with- 
out the  mass  can  be  obtained  from  Reference  27.  The  value  for  the  mass 
added  panel  damping  must  be  estimated  as  damping  data  for  panels  with  a 
heavy  center  mass  are  not  readily  available. 


3.3  Combined  Acoustic  and  Flight  Load  Panel  Stresses 


In  this  analysis , the  axially  applied  flight  loads  are  assumed  to  introduce 
a state  of  plane  stress  in  the  test  panel.  The  lateral  stress  from  the 
flight  loads  is,  therefore,  zero.  The  instantaneous  axial  stress,  oxF(t), 
which  is  assumed  to  be  uniform  over  the  width  of  the  panel  is,  therefore, 
related  to  the  instantaneous  axial  strain,  exF(t),  by 


°xF(t)  ' E<*F(t> 
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At  the  panel  edge,  the  axial  strain  from  the  simulated  acoustic  loading  is 
zero  (equation  15).  Consequently,  the  instantaneous  panel  edge  axial  and 
lateral  stresses,  axAF^)  3X1(1  0 respectively,  as  a result  of  the 

combined  shaker  simulated  acousticyand  flight  loads,  are  related  to  the  in- 
stantaneous axial  flight  load  stress,  and  the  instantaneous  lateral 

simulated  acoustic  load  stress,  ayA(t),  by 


°xAF^ 


qxF(t> 


U°yA(t) 


29 


ayAF^ 


o A(t) 

yAv 
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The  simui^cd  acoustic  loading  produces  a bending  stress  in  the  panel  skin 

5 maJ{imum  at  both  surfaces  and  zero  at  half  thickness.  The  surfaces 
of  the  skin  are,  consequently,  subjected  to  a biaxial  stress  at  at  1 L If 
interaction  effects  would  be  expected  to  initiate  from  either  surface  of  the 

ln  .identifyinS  interaction  effects  two  ratios  of  inplane  to  bend- 

noLl  tV+T  ? af  a milnimUm  Sh0uld  be  used‘  With  the  programmed  crack 
normal  to  the  inplane  load,  equations  29  and  30  indicate  that  a more  signi- 

f leant  effect  can  be  expected  by  changing  the  axial  flight  load  stress  While 
maintaining  a constant  simulated  acoustic  load  level. 

3 • Panel  Stress  Distribution 

In  the  previous  sections,  the  emphasis  is  placed  on  predicting  the  stress 
at  the  edge  of  the  panel  due  to  simulated  acoustic  loading.  As  long  as  the 
crack  propagates  along  the  panel  edge,  the  stress  at  the  panel  edge^ue  £0 
simulated  acoustic  loading  can  be  used  to  predict  the  stress  intensity  factor 
However,  if  the  crack  propagates  free,  the  edge  to  the  center  of  Se  pLelf 

“ ^\TC  r Pe^  stresses  at  the  edSe  which  diminish  to  zero  pLt  way 
in  and  then  peak  again  at  the  panel  center.  y 

ILtV  tbarefore>  necessary  to  define  the  stress  distribution  within  the 
panel.  The  curvature  measured  in  Reference  1 6 cannot  be  used  to  obtain  the 
stress  distribution  as  fixed-fixed  edge  conditions  were  simulated  in  SoJe 


In  the  U J ? attachment  for  the  panel  along  the  fastener  line. 

L ® d^ty\tb®  Peak  panel  edge  stress  location  is  influenced  by  the  panel 
edge,  especially  with  asymmetric  stiffness,  where  it  gravitates  towards  the 
hard  corner  of  the  stiffener  (Reference  23).  The  degree  of  edge  nSibmtv 

W y T accurately-  Since  empirical  data  indicates  sJmiS  7 

lateral  stress  magnitudes  for  the  edge  and  the  center  of  the  panel  a funda- 
mental modal  displacement  function  given  by  0.5  (l-  Cos  2 vy/b)  is^some- 
times  used  (Reference  19)  to  represent  the  lateral  panel  modal  deflection 
For  large  aspect  ratio  panels,  a simply  supported  axial  mode  shape  can  be* 
assumed.  The  above  mode  shapes,  when  used  in  the  Rayleigh  Ritz  analysis, 

1 produce  an  over- estimation  of  the  corresponding  resonant  frequency. 

diirihu+1  nnd^+i!aPeLCan  b<\ USed  t0  describe  both  the  strain  and  stress 
stribution  within  the  panel  in  conjunction  with  the  magnitudes  of  the 

strains  and  stresses  determined  by  the  analysis  in  previous  sections  or 
which  are  empirically  determined  from  measurements  at  the  panel  edge 

axial  and  lateral  panel  strains,  Tx(x,y)  and  fy^yj^esp^ely, 

at  point  x,  y,  are  given  in  terms  of  the  peak  lateral  panel  edge 
strain,  7 >"r  0 


ye- 


by 

(x,y)  = -0.25 


ye 


(» ) t1  ‘ 005  (t?)]  s1"  (it) 
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(x,y) 


COS 
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and  the  corresponding  stresses,  ax(x,y)  and  o (x,y),  respectively,  in 
terms  of  the  lateral  peak  panel  edge  stress,  dye,  hy 


(x,y)  = 


o (x,y)  = 0 

y ye 


-0-«(4-)2[l-coS(^)]  Sin  (I£)  + .cos  (S?)-.(a) 

- (¥) (?)  -0-25  u (?f  (¥)] sl"  (?) 
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4.  Simulation  of  Acoustic  Loading 

The  degree  of  acoustic  loading  simulation  possible  with  an  electromagnetic 
shaker,  apart  from  the  previously  derived  effects  on  panel  frequency  and 
stress,  can  be  determined  by  considering  the  generalized  load  or  force 
spectrum  given  by  equation  20  for  both  excitations.  The  increased  damping 
due  to  back  EMF  can  only  be  determined  experimentally.  Again  for  simpli- 
city, a simply  supported  panel  mode  shape  is  assumed  in  the  subsequent 
analysis. 

4.1  Acoustic  Excitation 


Two  types  of  acoustic  excitation  are  usually  considered,  namely;  normal 
incidence  and  grazing  incidence.  For  the  fundamental  panel  mode  with 
normal  incidence  noise  and  even  with  grazing  incidence  noise  when  the 
acoustic  wavelength  is  considerably  greater  than  the  panel  length,  the 
pressure  field  is  fully  correlated  over  the  surface  of  the  panel.  The 
generalized  load  spectrum,  equation  17,  reduces  to  the  separate  integral 
of  the  mode  shapes  over  the  surface  of  the  panel  multiplied  by  a constant 
noise  spectrum  level,  Gq.  For  the  fundamental  mode  of  a simply  supported 
panel,  the  fully  correlated  generalized  load  spectrum  is  reduced  to 


GFa 


16a2  b2 
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With  grazing  incidence  excitation  for  higher  order  panel  modes  and  for  the 
fundamental  mode  when  the  acoustic  wave  length  approaches  the  panel  length, 
the  expression  for  the  generalized  load  spectrum  (Reference  2 6)  is 


8a2  b2 

f1  ±cos(?) 

gfa  - | 

Pi 

2 

2 2 

,2  - M 

\cm/  -1 

n m | 
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where  + is  used  with  odd  m,  n 

- is  used  with  even  m and  odd  n 

If  the  longest  panel  side  is  parallel  with  the  acoustic  propagation  path, 
it  is  possible  to  excite  the  (2,l)  mode,  especially  for  high  aspect  ratio 
panels  for  which  the  (2,l)  mode  frequency  can  approach  the  frequency  of 
the  fundamental  mode.  The  (2,l)  mode,  in  addition,  has  a very  low  damping 
coefficient  due  to  acoustic  radiation  cancellation  effects.  The  (2,1)  mode 
can,  however,  be  eliminated  by  rotating  the  panel  with  the  long  side  normal 
to  the  propagation  path.  The  (l,2)  mode  that  could  now  be  excited  would 
have  a considerably  higher  frequency  than  the  fundamental  mode.  Also,  the 
(3,l)  mode,  after  the  panel  is  rotated,  will  have  a small  generalized  load 
spectrum  in  comparison  to  the  fundamental  mode. 

k.2  Shaker  Excitation 

Even  modes  of  the  panel  cannot  be  excited  by  the  shaker  when  it  is  attached 
to  the  center  of  the  panel.  The  shaker  will,  however,  excite  all  the  odd 
modes  of  the  panel  equally  when  a constant  force  spectrum  level  is  used. 
With  acoustic  excitation,  a smaller  generalized  load  spectrum  (equation  36) 
is  generally  obtained  for  the  higher  order  panel  modes.  Unlike  the  acous- 
tic excitation,  the  shaker  is  a point  drive  system  acting  at  the  center  of 
the  panel  where  the  mode  shape  is  unity.  Thus,  the  shaker  generalized 
force  spectrum  reduces  to 


GFS  = FS 
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where  Fg  is  the  constant  shaker  force  spectrum  level. 

^ • 3 Ratio  of  Shaker  to  Acoustic  Force 

The  generalized  load  spectrum  for  acoustic  excitation  can  be  written  in  terms 
of  the  force  spectrum,  F^,  acting  on  the  panel  and  given  by 

Fa  = a2b2  00  38 

The  ratio  of  the  shaker  generalized  force  spectrum  to  acoustic  generalized 
load  spectrum  becomes 


2!s  ’LZs 

GF.  16F.  39 

A A 

Thus,  the  shaker  provides  greater  excitation  for  the  same  force  spectrum 
level  than  does  the  acoustic  excitation. 
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4.4  Shaker  Damping 


If  the  ratio  of  the  shaker  coil  weight  to  the  panel  weight  is  large,  it  is 
to  be  expected  that  the  damping  from  the  shaker  due  to  back  EMF  effects 
would  also  be  large.  A damping  coefficient  in  the  order  of  0.4  could  be 
possible  under  these  circumstances.  The  effect  of  damping  on  the  normalized 
response  spectrum  of  a single  degree-of-freedom  system  is  illustrated  in 
Figure  3»  The  normalized  single-degree-of-freedom  response  spectrum  ,H(— ) , 
obtained  by  multiplying  equation  18  by  Mr2  wr4,  is  given  by  \rl 


With  a damping  coefficient  of  0.4,  the  response  peak  is  shifted  towards  a 
lower  frequency  while  the  amplification  near  resonance  is  almost  constant 
and  close  to  unity.  Under  these  circumstances  it  would  be  almost  impossible 
to  experimentally  detect  the  resonant  frequency  and  to  obtain  significant 
panel  response  with  broad  band  random  excitation. 

4 . 5 Band  Limited  and  Forced  Excitation 

Limiting  the  bandwidth  of  the  excitation  to  increase  the  excitation  spectrum 
level  has  been  widely  used  in  sonic  fatigue  testing  to  obtain  a greater 
stress  level  in  the  test  specimen.  The  effect  of  limiting  the  excitation 
bandwidth  on  the  mean  square  response  of  a single  degree-of-freedom  system 
has  been  determined  from  the  theory  in  Reference  28  and  is  presented  in 
Figure  4,  normalized  by  the  mean  square  response  from  white  noise  excita- 
tion. In  general,  the  excitation  bandwidth  has  to  be  at  least  3»3  times  the 
3dB  point  bandwidth  for  a loss  of  10  percent  or  less  in  the  rms  response 
level.  This  method  is  effective  for  structures  with  relatively  low  damping 
coefficients  and  correspondingly  high  amplification  factors.  In  the  case 
of  high  damping  coefficients  producing  a broad  panel  response  and  near 
unity  amplification,  high  stress  levels  can  only  be  obtained  by  artificially 
increasing  the  amplification  factor  through  the  use  of  very  narrow  band 
excitation.  Narrow  band  filters  have  a shape  similar  to  a single  degree-of- 
freedom  system  over  the  significant  portion  of  the  spectrum  and,  therefore, 
can  be  used  to  simulate  the  response  of  the  single  degree-of-freedom  system. 
The  filter  bandwidth  should  be  as  close  as  possible  to  the  desired  single 
degree-of-freedom  response  3dB  bandwidth. 

In  acoustic  fatigue  tests  and  crack  propagation  tests,  the  crack  growth  and 
sonic  fatigue  life  are  a function  of  the  rms  panel  stress  (References  l6 
and  19) . A Rayleigh  distribution  for  the  stress  peaks  is  generally  obtained 
in  the  test  panel.  A narrow  band  forcing  function  would  produce  similar 
effects  in  the  test  panel.  A Rayleigh  distribution  of  stress  peaks  is 
obtained  with  an  ideal  rectangular  narrow  band  filter,  indicating  that  the 
filter  shape  is  not  a significant  parameter. 
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FIGURE  3.  NORMALIZED  SINGLE  DEGREE  OF  FREEDOM 
FREQUENCY  RESPONSE  FUNCTION 
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The  use  of  a filter  simulated  response  function  will  make  it  difficult  to 
simulate  the  variation  of  panel  response  frequency  with  the  inplane  load 
normally  encountered  in  structures  with  low  damping  when  excited  by  broad 
band  noise.  With  a highly  damped  response,  a reasonable  simulation  can  be. 
achieved  with  a constant  excitation  frequency.  Some  gain  in  response  ampli- 
tude can  be  obtained  by  employing  a lower  excitation  frequency  than  the 
panel  fundamental  frequency  (zero  damping). 

5.  Detailed  Design  of  Test  Specimen 

The  final  design  of  the  test  specimen  illustrated  in  Figure  ^ is  a compro- 
mise between  the  requirements  of  the  compressive  load  carrying  capability  of 
the  panel,  the  interaction  of  the  inplane  loads  with  the  centrally  attached 
mass  panel  modal  frequencies,  the  panel  stresses  due  to  shaker  and  acoustic 
excitations,  and  the  capacity  of  the  MTS  fatigue  machines  for  accurately 
controlling  the  tape  load  reproduction.  This  last  condition  imposed  a mini- 
mal cross  sectional  area  requirement  on  the  test  specimen. 

5.1  Test  Specimen  Material  and  Thickness 

7075-T6 • aluminum  alloy  is  used  as  the  material  for  the  panel,  skin,  doublers 
and  the  frames.  The  need  for  a small  test  specimen  cross  sectional  area  and 
adequate  depth  for  the  countersunk  fastener  hole  resulted  in  the  employment 
of  a O.063  inch  panel  thickness.  The  extruded  channel  section  stiffener  has 
a thickness  of  0.07  inches,  a height  of  3 inches  and  a flange  width  of  O.85 
inches.  Good  sonic  fatigue  practice  requires  the  stiffener  to  be  one  gage 
thicker  than  the  panel,  while  the  frame  height  was  selected  to  ensure  a 
sufficiently  high  frame  bending  frequency  relative  to  the  panel  fundamental 
frequency,  to  preclude  coupling  effects. 

5.2  Selection  of  Panel  Size 

Preliminary  analysis  indicates  that  a panel  compression  load  stress  and  a 
critical  panel  buckling  stress  of  7500  psi  and  15,000  psi,  respectively, 
would  be  realizable  with  the  test  specimen  while  meeting  the  other  design 
requirements.  The  7500  psi  difference  is  the  safety  margin  used  to  pre- 
clude panel  buckling, due  to  the  uncertainty  in  the  panel  edge  stiffness. 


For  a given  panel  thickness,  the  magnitude  of  the  critical  panel  buckling 
stress  essentially  determines  the  panel  width.  The  critical  panel  mode  is, 
however,  dependent  on  the  panel  length.  Ih  panels  with  large  aspect  ratios, 
the  higher  order  panel  modes  produce  the  lowest  critical  buckling  stresses. 
A large  panel  aspect  ratio  is,  however,  desirable  to  obtain  a near  constant 
panel  edge  stress  distribution  near  the  center  of  the  longest  panel  edge  in 
the  fundamental  mode.  Unfortunately,  higher  order  odd  panel  modes  would 
also  be  excited  by  band  limited  random  excitation  needed  to  accomodate  the 
theoretically  predicted  panel  fundamental  frequency  variation  due  to  the 
inplane  flight  loads,  at  large  panel  aspect  ratios. 
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NOTE:  ALL  DIMENSIONS  ARE  IN  INCHES 
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FIGURE  5.  TEST  PANEL  DETAILS 


A panel  width  of  4.5  inches  and  length  of  10  inches  are  selected  on  a basis 
of  calculations  employing  an  interative  type  procedure  which  is  schemati- 
cally illustrated  in  Figure  6 . The  calculated  panel  frequency  variation 
with  inplane  panel  stress  and  aspect  ratio  are  illustrated  in  Figures  7 
and  8,  respectively.  Equation  11,  which  is  used  to  determine  the  panel 
frequency,  may  not  be  accurate  for  the  higher  order  axial  panel  modes  due 
to  the  greater  influence  of  the  attached  panel  center  mass  on  the  higher 
order  axial  mode  shapes. 

The  analysis,  however,  indicates  the  possible  occurrence  of  a 3,1  mode  in 
addition  to  the  fundamental  mode  at  a compressive  stress  of  7500  psi.  The 
large  shaker  coil  weight  of  3.09  pounds  is  responsible  for  the  relatively 
low  but  acceptable  panel  fundamental  frequency. 

Panel  edge  crack  growth  data  from  random  acoustic  excitation  (Reference  l6) 
indicates  that  a stress  level  in  the  region  of  3300  psi  rms  would  be  neces- 
sary at  the  panel  edge  to  produce  a crack  growth  rate  commensurate  with  the 
available  test  time.  Pretest  analysis  of  the  shaker  attached  panel  response 
indicates  that  with  the  level  of  force  available  in  the  acoustic  test  facil- 
ity, it  is  possible  to  reach  3300  psi  rms  stress  level  with  a damping  coeffi- 
cient below  2 percent.  However,  spectrum  shaping  is  necessary  to  obtain  an 
increased  spectrum  level  around  the  panel  resonance. 

No  pretest  analysis  was  conducted  for  the  narrow  band  filtered  force  simula- 
tion of  the  single  degree-of -freedom  system  response  employed  in  the  test 
program  to  overcome  the  extremely  high  shaker  damping.  The  magnitude  of 
the  shaker  damping  could  only  be  established  during  the  calibration  tests, 
thus  the  simulation  had  to  be  developed  experimentally. 

5.3  Design  Detail 

The  panel  ends  are  designed  to  transfer  the  inplane  loads  from  the  end  grips 
by  friction  with  the  least  possible  stress  concentration  effect.  The  change 
in  panel  width  employing  large  radii,  together  with  the  staggered  doublers 
are  designed  to  produce  a uniform  distribution  of  the  axial  panel  stress 
across  the  panel  width. 

To  ensure  that  the  crack  propagation  is  confined  to  the  panel  between  the 
fastener  lines,  a free  panel  edge  extending  beyond  the  frame  edge  is  em- 
ployed. Analyses  were  carried  out  for  the  panel  free  edge  and  the  stiffen- 
er frame  flange  buckling,  the  inter  rivet  buckling  and  the  long  column 
frame  buckling  up  to  the  critical  panel  buckling  stress.  Panel  fastener 
spacing  was  also  checked  for  sonic  fatigue  requirements  while  fastener  edge 
distances,  compatible  with  good  fatigue  design,  are  employed  throughout 
the  test  specimen.  Eight  3/l6  inch  diameter  countersunk  fasteners  are 
employed  around  the  center  of  the  long  panel  sides  to  represent  actual  air- 
frame construction.  The  shallow  head  countersunk  fasteners  are  selected 
to  avoid  knife  edges  in  the  fastener  holes. 


21 


FIGURE  6.  ITERATIVE  PANEL  DESIGN  PROCEDURE 


22 


MODE  1,1 

MODE  2,1 

MODE  3,1 


PANEL  DIMENSIONS 

a = 10" 

b = 4.5" 
h = .063" 


AXIAL  INPLANE  STRESS  - KSI 


FIGURE  7.  EFFECT  OF  AXIAL  INPLANE  STRESS  ON  PANEL  FREQUENCY 
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FIGURE  8.  EFFECT  OF  ASPECT  RATIO  ON  FREQUENCY  OF  PANEL  SUBJECTED  TO  COMPRESSION 

STRESS  OF  7500  PS I. 


The  B&K  force  gage  is  bolted  as  close  to  the  panel  surface  as  practical 
to  measure  the  true  force,  free  of  flexibility  effects  (Figure  9).  The 
panel  center  is  reinforced  by  bonding  a 0.032  inch  thick  circular  doubler 
(0.75  inch  diameter)  to  each  side  of  the  panel  per  Lockheed  specification 
LCP79- 206l-Type  20.  A washer  is  inserted  between  the  doubler  and  the 
force  gage  to  provide  a more  gradual  stress  transfer  to  the  panel.  The 
shaker  drive  rod  is  attached  to  the  other  side  of  the  force  gage  through  a 
machined  adapter.  Analysis  of  the  shaker  rod  bending  frequency  was  con- 
ducted to  assist  in  reducing  the  probability  of  resonant  coupling  with  the 
panel  fundamental  mode. 

5.4  End  Grip  Design 

Existing  MTS  fatigue  machine  end  grips  are  modified  by  machining  the  outer 
end  grip  surfaces  parallel  to  the  panel  axis.  The  upper  attachments 
(Figure  10)  are  machined  from  one  inch  thick  steel  plate  to  fit  on  the 
outer  surface  of  the  end  grips  and  pick  up  the  frame  loads  through  a full 
length  doubler  on  each  frame  flange.  The  lower  attachments  are  machined 
from  l/4  inch  thick  steel  plate  to  the  same  shape  as  the  panel  ends  and  fit 
over  the  panel  ends  to  provide  a friction  type  grip  with  a gradual  load 
transfer  into  the  face  sheet.  The  test  panel  ends  fit  over  the  other  outer 
surface  of  the  end  grips.  This  arrangement  is  dictated  by  the  need  for  the 
load  axis  to  pass  through  the  centroid  of  the  test  specimen  cross  section 
to  eliminate  frame  bending  and  the  associated  panel  stress. 

5 • 5 Fatigue  Crack  Initiation 

A single  simulated  crack  is  employed  normal  to  the  panel  fastener  line  in 
one  of  the  countersunk  fastener  holes  (Figure  11 ) nearest  to  the  center 
of  the  panel,  in  all  test  panels  except  the  calibration  panel.  The  simu- 
lated crack  is  made,  initially,  by  cutting  a slot  with  a jeweler's  saw 
into  the  face  sheet  at  the  selected  fastener  hole,  prior  to  final  assembly. 
After  assembly,  with  the  test  specimen  mounted  in  the  fatigue  machine,  the 
test  specimen  is  fatigue  loaded  at  a reduced  level  relative  to  the  test 
load  until  a 0.01  inch  long  crack  is  obtained  at  the  end  of  the  slot.  An 
additional  simulated  crack  is  made  parallel  to  the  fastener  line  (Figure  12) 
in  the  panel  center  fastener  hole  adjacent  to  the  above  crack,  in  two  of 
the  test  panels  subjected  only  to  shaker  simulated  acoustic  loading.  The 
additional  crack  location  is  typical  of  a sonic  fatigue  induced  crack. 

Crack  initiation  in  this  instance  is  by  the  application  of  shaker  simulated 
random  loading. 
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SHAKER  DRIVE  ROD 


FIGURE  9.  FORCE  GAGE 
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FIGURE  10.  PANEL  END  GRIP  ASSEMBLY 


SECTION  AA' 


FIGURE  11.  PRECRACK  FOR  COMBINED  LOADS 


SECTION  BB' 


FIGURE  12.  PRECRACK  FOR  SHAKER  SIMULATED 
ACOUSTIC  LOADING 
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Ill  LOAD  HISTORIES  AND  SPECTRA 


1.  Introduction 


The  discussion  in  this  section  describes  the  development,  from  existing  data, 
representative  flight  load  and  acoustic  test  spectra  and,  thereafter,  the 
corresponding  time  histories  which  are  recorded  onto  magnetic  tares.  The 
flight  load  tape  typically  contains  pseudo  random  flight-by- flight  load  time 
histories  of  representative  runway,  landing  impact,  gust  and  maneuver  load 
spectra.  The  acoustic  tape  contains  typical  takeoff  and  accelerated  flight 
jet  noise  time  histories.  These  tapes  will  be  used  to  control  the  load 
applied  to  the  test  panels  by  both  the  MTS  fatigue  machines  and  the  electro- 
magnetic shakers. 

The  flight  load  spectra  along  with  the  time  histories,  in  general,  have  to 
be  modified  to  conform  to  the  requirements  imposed  by  the  specimen  design, 
test  facility  capability  and  the  load  tape  duration.  Truncation  of  the  low 
level  portion  of  the  load  spectra,  for  example,  is  carried  out  due  to  the 
limitation  on  number  of  flight  cycles  that  can  be  accommodated  within  a 
reasonable  flight  duration.  This  truncation  permits  the  application  of  high 
amplitude  cyclic  loads  at  lower  frequency  and  with  high  accuracy. 

The  real  time  simulation  of  the  ground  acoustic  loading,  in  this  program, 
amounts  to  a duration  of  approximately  18  seconds  of  the  total  acoustic  loading 
of  each  flight.  The  duration  of  the  flight  portion  is  determined  by  the 
flight  load  requirements.  The  lower  level,  but  very  much  longer  duration 
of  the  flight  acoustic  load,  can  introduce  a significant  effect  on  crack 
growth  which  is  spread  over  a longer  period  of  time.  To  compensate  for  this 
situation,  the  level  of  the  flight  acoustic  loading  can  be  increased  in  a 
trade-off  against  time,  to  conform  to  the  flight  duration  imposed  by  the  flight 
spectral  requirements . 

The  analysis  in  Section  II  indicates  that  at  least  two  ratios  of  flight  loads 
to  acoustic  loads  are  desirable  to  identify  trends.  The  propulsion  config- 
urations in  advanced  short  takeoff  and  landing  (STOL)  transport  type 
airplanes  indicate  that  both  upper  and  lower  wing  flight  loads  in  combination 
with  acoustic  loads  are  of  interest.  The  upper  and  lower  wing  spectra  are 
generally  of  different  magnitudes.  Thus  the  employment  of  the  upper  and 
lower  wing  load  spectra  with  a constant  level  acoustic  excitation  will 
satisfy  both  of  the  above  needs . 

2.  Flight-by-Flight  Loads 


The  flight-by-flight  load  spectra  are  generally  based  on  mission  profiles 
and  operation  load  spectra  for  a specific  aircraft  category.  Such  spectra 
indicate  the  relative  frequency  but  not  necessarily  the  sequence  with  which 
specific  loading  levels  are  reached  or  exceeded  during  a specified  period  of 
time.  The  mission  profile  defines  the  number  and  duration  of  each  mission 
type  and  the  number  of  landings  for  the  life  of  the  aircraft.  This  data  estab- 
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lishes  the  proper  mix  ratio  to  he  used  in  the  load  tapes.  The  individual 

alternating  load  cycles  within  each  mission  is  selected  randomly  from  the 
pool  of  appropriate  load  spectra  corresponding  to  the  specific  mission. 

The  accumulated  number  of  load  cycles  for  each  mission  category  is  based  on 
the  load  spectra  and  the  duration  of  both  the  acoustic  and  flight  loadings 
for  each  flight. 

2.1  Basic  Mission  and  Spectral  Data 

The  basic  mission  profiles  and  the  flight-by-flight  load  spectral  data 
employed  in  this  program  have  been  obtained  from  References  29  and  30  for  a 
STOL  type  airplane.  The  simplified  three  mission  profiles,  for  the  STOL 
airplane,  are  summarized  in  Table  1 and  schematically  illustrated  in  Figure 
13.  The  corresponding  spectra  for  these  missions  are  summarized  in  Tables 
2 through  5.  Table  2 represents  the  taxi,  takeoff  and  landing  roll  spectra. 
The  maneuver  and  gust  load  spectra  (Table  3(a))  and  the  corresponding 
turbulance  parameters  (Table  3(b))  are  based  on  MIL-A-8866A,  Tables  VII  and 
VIII.  Table  4 contains  the  gust -plus -maneuver  spectrum  for  the  low  level 
resupply  mission  and  Table  5 contains  the  landing  impact  spectrum. 

As  a further  simplification,  the  three  mission  catagories  (Table  l)  are 
reduced  to  two,  namely,  the  typical  and  the  low  level  resupply  missions, 
by  combining  the  basic  training  missions  with  the  typical  missions.  The 
duration  of  each  flight  within  the  two  missions  has  been  fixed  at  0.5 
hours  and  0.9  hours  for  the  low  level  resupply  and  the  typical  mission, 
respectively.  Since  the  number  of  typical  missions  is  approximately  three 
times  that  of  the  low  level  resupply  missions,  a total  number  of  40,250 
takeoffs  and  landings  are  calculated  for  the  30,000  hour  life  of  the 
airplane  using  the  above  flight  durations.  Touch-and-go  landings  have 
been  omitted. 

Simulated  flight  duration  of  approximately  28  seconds  takes  into  account 
an  18  second  duration  for  the  takeoff  accoustic  loading  and  approximately 
a 10  second  duration  for  the  flight  load.  The  number  of  flights  that  can 
be  accomodated  on  a 90  minute  tape  is,  therefore,  around  185  which  also  allows 
time  for  the  initial  calibration  signal.  This  permits  a selection  of  74 
missions  of  which  18  missions  are  considered  to  be  in  the  low  level  resupply 
mission  category  and  the  remaining  56  missions  in  the  typical  mission 
category.  Table  1 and  Figure  13  indicate  that  each  low  level  resupply  mission 
has  four  flight  segments, while  a typical  mission  has  two  flight  segments. 
Hence,  the  actual  number  of  flights  that  can  be  accomodated  by  a 90  minute 
duration  load  tape  is  184. 

The  actual  duration  of  the  ground  and  flight  section  in  each  flight  does 
vary  for  the  upper  and  lower  wing  surface  loading  and  differs  slightly 
from  the  above  values  in  the  final  tape.  The  intention,  however,  is  to 
maintain  the  same  values  for  these  flight  sections  in  both  upper  and  lower 
wing  surface  load  tapes. 
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FLT.  NO. 


FLT.  NO. 


FLT.  NO. 


FLT.  NO. 


FLT.  NO. 


FLT.  NO. 


1-1A 


1-1B 


1-1C 


1-2A 


1-2B 


2-1A 


2-2A 


J 


500  FT  .\\J  500  FT.  \ / 500  FT.  \ / 500  FT. 

STOL  STOL  STOL  CTOL 

3-2 


500  FT 


\J  500  FT  500  FT.y^ 


500  FT.) 


CTOL  CTOL  CTOL  CTOL 

NOTE:  Figure  Taken  From  Reference  30. 


FIGURE  13.  PABST  AIRCRAFT  UTILIZATION  (SIMPLIFIED) 
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TABLE  2.  TAXI,  TAKEOFF,  AND  LANDING  ROLL 
SPECTRUM  (PREPARED  AND  SEMI- 
PREPARED  RUNWAYS) 


± Ag 

CUMULATIVE  OCCURRENCES/lOOO  TAXI'S 

Prepared  Runways 

Semi-Prepared  Runways 

Unprepared  Runways 

.1 

194,094 

280,000 

280,000 

.2 

29,094 

100,000 

210,000 

.3 

2,094 

30,000 

130,000 

.4 

94.155 

5,000 

68,000 

.5 

4.155 

750 

22,000 

.6 

0.155 

100 

5,300 

.7 

0.005 

15 

900 

.8 

2 

130 

.9 

0.3 

20 

1.0 

0.05 

3 

l.i 

0.5 

NOTE:  Table  Taken  From  Reference  30,  Attachment  3. 
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TABLE  3(a).  MANEUVER  LOAD  SPECTRA  FOR 
TRANSPORT  AIRCRAFT 


Spectrum  Equation 

1.  SN(Ag)  = Nq1  eAs/bl  + Nq2  eAS/b2 


2.  EN(Ag)  = Nq1  e 


~^S2/2b^2 


+ Nq2  e 


-Ag2/2b22 


Mission 

Type 

EQ. 

# 

N *i 
ol 

OJ 

o 

& 

bl 

b2 

CAssault 

* 

Climb 

+ 

1 

2.7  x 105 

2.5  x 102 

.0492 

.1646 

+ 

1 

9.0  x io4 

50.0 

.0869 

.2877 

Cruise 

+ 

1 

, 4 

4.o  x io 

103 

.0523 

.0986 

+ 

1 

4 

2.1  x 10 

62.0 

.0921 

.2621 

Descent 

+ 

1 

1.9  X 105 

0 

.0543 

0 

+ 

1 

1.7  X 10b 

4.5  x 102 

.0942 

.2311 

C 

Transp. 
(Logistics ' 

Climb 

+ 

2 

9.0  x 103 

103 

.1200 

.1614 

+ 

1 

2.6  x 10b 

54 

.06157 

.2942 

Cruise 

+ 

2 

4.52x  io2 

95 

.1452 

• 

H 

00 

OO 

VJ1 

+ 

1 

9.0  x io1* 

75 

.o4o8 

.1687 

Descent 

+ 

2 

8.83x  103 

2.2  x 102 

.1336 

.1827 

+ 

1 

3.5  x 10^ 

4.0  x io2 

.05509 

.19054 

Cm 

Transp. 

(Training) 

** 

Climb 

+ 

1 

1.6  x 105 

0 

.0598 

0 

+ 

5.2  x 10b 

2.7  x 103 

.0835 

.1820 

Cruise 

+ 

IT~ 

4.0  x 10 

0 

.0566 

0 

+ 

4.6  x 10b 

4o 

on 

00 

0 

• 

.2966 

Descent 

+ 

4 

5.5  x 10 

9.5  x 103 

.0476 

.0841 

+ 

1 

5.4  x 105 

1.6  x 10^* 

.0670 

.1922 

* MIL-A-8866A 
**  MIL-A-8866A 


Table  VIII  Assault  Spectrum  for  Transport  Airport 

Table  VII  Transport  and  Training  Spectra  For 

Transport  Aircraft 


NOTE:  Table  Obtained  From  Reference  30,  Attachment  3. 
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TABLE  3(b).  TURBULENCE  PARAMETER 
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TABLE  4.  LOW  LEVEL  RE-SUPPLY  MISSION 

GUST  PLUS  MANEUVER  LOAD  SPECTRIN 

Alt.  - 1000  Ft. 

V = 300  KEAS 


± Ag  c.g. 

Cumulative  Occurrence 
Per  Flight  Hour 

.1 

1,300 

.2 

420 

.3 

160 

.4 

52 

.5 

17 

.6 

4 

.7 

2.8 

.8 

0.95 

NOTE:  Table  Taken  From  Reference  30,  Attachment  3» 
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TABLE  5.  LANDING  IMPACT  INCREMENTAL  LOAD 
SPECTRUM 


Ag 

CTOL  Per 
1,000  Landings 

STOL  Per  1,000  Landings 

S 

G.W. < 153,OOOff 

3 

G.W.  > 153, 000^ 

.1 

530 

.2 

50 

900 

800 

.3 

9 

600 

290 

.4 

4.3 

450 

160 

.5 

2.2 

320 

94 

.6 

1.3 

250 

60 

.7 

0.65 

210 

39-1 

.8 

163 

25.02 

• 9 

133 

14.04 

1.0 

110 

10.99 

1.2 

70 

4.58 

1.4 

50 

1.6 

37 

1.8 

26 

2.0 

20 

2.2 

15.5 

2.4 

11.3 

2.6 

7.4 

2.8 

4.8 

3.0 

3.2 

3.2 

2.2 

3.4 

NOTE:  Table  Taken  From  Reference  30,  Attachment  3. 
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2.2  Lover  Wing  Surface  Load  History 

Except  for  portions  of  the  takeoff  and  landing  runs,  the  flight  conditions 
will  impose  tension  loads  on  the  lower  wing  surface.  The  flight-by- flight 
load  history  is  generated  using  the  following  described  techniques.  A 
total  of  50  load  cycles  can  be  accommodated  within  each  flight  segment  when 
using  the  preselected  flight  duration  of  10  seconds  and  a nominal  loading 
frequency  of  5 cps.  There  are  72  flight  segments  for  the  18  low  level 
resupply  missions.  The  total  number  of  load  cycles  is  then  3600  cycles. 

The  truncated  load  spectrum  is  shown  in  Figure  14  together  with  a spectrum 
based  on  Table  4 with  18  missions  (36  flight  hours).  The  number  of  cycles 
at  selected  load  levels  are  also  shown  in  Table  6.  The  load  spectra  for 
the  typical  mission  and  the  low  level  resupply  mission  are  presented 
seperately.  The  high  load  level  truncation  is  used  to  eliminate  the 
infrequent  high  load.  The  low  load  level  truncation  is  based  on  the 
limit  of  the  total  number  of  3600  cycles  allotted  to  the  18  low  level 
resupply  missions.  The  two  high  and  low  load  level  truncation  justifications 
apply  equally  for  the  other  spectra. 

Similar  truncation  is  carried  out  for  the  spectrum  of  the  56  typical  missions. 
The  typical  mission  has  two  flight  segments.  Based  on  50  cycles  per  flight 
segment,  the  total  number  of  load  cycles  is  56OO  cycles  which  is  more  than 
enough  to  accommodate  the  whole  spectrum  as  given  in  Table  3.  However,  the 
low  load  level  truncation  atAg  = 0.2  is  selected  and  2600  cycles  of  Ag  = 0.0 
are  inserted  to  replace  the  surplus  cycles.  The  high  load  level  truncation 
is  atAg  = 0.7  which  is  the  same  as  for  the  low  level  resupply  missions. 

The  resulting  truncated  spectrum  is  shown  in  Figure  14  and  Table  7. 

The  taxi,  takeoff,  and  landing  roll  spectrum  is  reduced  to  a constant 
sustained  load  covering  the  duration  of  the  entire  ground  load  segment  of  the 
derived  spectrum.  This  simplification  can  be  justified  by  the  reasoning 
that  compressive  cyclic  loads  about  a compressive  mean  load  do  not  have 
significant  effect  on  the  crack  growth  behavior.  The  value  of  the  sustained 
load  factor  is  selected  to  be  0.5g  based  on  the  critical  buckling  load  of 
the  test  panel. 

From  Table  5>  it  can  be  seen  that  the  once  per  landing  incremental  load 
factor  (for  landing  impact)  is  approximately  0.6  to  0.7g.  This  value  is 
similar  to  the  sustained  load  level  considered  for  the  ground  load  segment. 
Hence,  the  landing  impact  spectrum  is  not  considered  in  the  generation  of 
the  load  history. 

The  flight -by-flight  load  history  is  generated  from  the  modified  load 
spectra  using  a computerized  randomization  procedure.  From  the  total 
of  l8  low  level  resupply  missions  and  56  typical  missions,  7 mission  types 
composed  of  4 low  level  and  3 typical  mission  types  are  selected  for 
randomized  repetition.  The  sequences  of  the  individual  load  cycles  within 
each  of  the  four  flight  segments  in  a low  level  resupply  mission  and  the 
three  flight  segments  in  the  typical  mission  are  randomized.  The  number  of 
repetitions  for  each  of  the  7 mission  types  and  the  number  of  individual 
load  cycles  within  each  mission  type  are  selected  to  match  as  close  as 
possible  the  flight  load  spectra  in  Figure  14.  Thus  in  the  low  level 
missions,  Mission  Types  1 and  2 are  repeated  five  times  each  and  Mission 
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FIGURE  14.  TRUNCATION  OF  LOAD  SPECTRA 


TABLE  6.  THE  TRUNCATED  LOW  LEVEL  RE-SUPPLY 
MISSION  SPECTRUM  FOR  l8  MISSIONS 


Ag 

No.  of  Cycles 

0.74 

72 

0.66 

78 

0.58 

150 

0.52 

300 

0.48 

400 

0.44 

800 

0.36 

1800 

TABLE  7.  THE  TRUNCATED  TYPICAL  MISSION 
SPECTRUM  FOR  56  MISSIONS 


Ag 

No.  of  Cycles 

0.70 

7 

0.64 

8 

0.56 

20 

0.50 

35 

0.44 

50 

0.36 

310 

0.26 

950 

0.16 

2100 

40 


Types  3 and  4 are  repeated  4 times  each*  In  the  typical  mission.  Mission 
Type  1 is  repeated  40  times  and  Mission  Types  2 and  3 are  repeated  4 times 
each.  The  number  of  mission  types  and  the  number  of  repetions  are  selected 
using  as  a guideline  the  number  of  load  cycles  at  the  high  truncation  load 
level.  The  loads  with  a frequency  of  occurence  less  than  once-per-flight 
are  generally  assigned  to  one  mission  type,  making  it  mors  severe. 

The  conversion  factor  from  lg  flight  load  to  stress  is  7400  psi  as  shown 
in  Reference  29.  The  ground  load  and  landing  impact  is  simulated  using  a 
sustained  constant  load  factor  of  -0.5g.  The  conversion  factor  from  lg 
ground  load  to  stress  is  5000  psi  (Reference  29). 

2.3  Upper  Wing  Surface  Load  History 

A composite  flight-by- flight  load  history  can  be  generated  also  for  the 
upper  wing  surface  following  a similar  procedure  to  that  described  above 
for  the  lower  wing  surface.  However,  significant  quantitative  differences 
exist  between  the  two  composite  load  histories  primarily  because  on  the 
upper  wing  surface  the  tension  cyclic  loads  are  induced  by  the  ground  loads 
and  the  compressive  cyclic  loads  are  induced  by  the  flight  loads.  This 
consideration  leads  to  the  selection  of  an  18  second  ground  load  segment 
and  a 10  second  flight  load  segment.  The  selection  of  a 10  second  flight 
duration  instead  of  a shorter  time  span  is  required  to  maintain  the  same 
flight  acoustic  load  duration  as  employed  in  the  lower  wing  surface  load 
history. 

Based  on  the  18  second  duration  on  the  ground  and  a nominal  loading  frequency 
of  5 cps,  the  total  number  of  cycles  obtained  for  the  184  ground  load  segments 
is  16,560  cycles.  The  taxi,  takeoff  and  landing  roll  spectrum,  shown  in  Table 
2,  is  modified  accordingly.  The  high  load  level  truncation  is  selected  to  be 
0.8g.  The  low  load  level  truncation  is  limited  by  the  total  number  (16,560) 
of  cycles  allowed  for  ground  load  segments. 

One  additional  step  was  included  in  constructing  the  taxi,  takeoff,  and 
landing  roll  spectrum  for  the  basic  set  of  72  missions  (a  total  of  184 
flight  segments)  to  reflect  a proper  mix  of  the  three  runway  configurations 
which  are:  (1)  unprepared  runway,  (2)  semi-prepared  runway,  and  (3) 

prepared  runway.  A runway  configuration  mix  for  the  three  runway  config- 
urations is  obtained  using  the  mission  profile  data  in  Table  1.  The  resulting 
runway  configuration  mix  includes  42  percent  of  runway  (l),  24  percent  of 
runway  (2),  and  34. percent  of  runway  (3).  Based  on  the  above  mix  of  runway 
configurations,  the  milder  spectrum  shown  in  Figure  14  and  in  Table  8, 
has  been  constructed.  The  tension  loads  induced  by  the  landing  impacts  are 
neglected. 

A similar  randomization  procedure  to  that  employed  for  the  lower 
wing  surface  is  used  to  generate  the  flight-by-flight  ground  load  time 
history.  Based  on  the  load  spectra  in  Table  8,  eight  basic  ground  load 
segments  are  selected,  each  repeated  23  times  in  a randomized  sequence. 

The  individual  load  cycles  within  each  segment  are  also  randomized. 


4l 


The  lg  stress  conversion  factor  is  equal  to  5000  issi.  The  somewhat 
arbitrary  sustained  load  factor  of  -lg  is  selected  for  the  flight  segments. 
This  value  corresponds  to  a compression  stress  of  5000  psi  which  is  just  below 
the  buckling  limit  calculated  for  the  test  panels. 

TABLE  8.  THE  TRUNCATED  TAXI,  TAKEOFF,  AND  LANDING 
SFECTRUM  FOR  72  MISSIONS 


Incremental 
Load  Factor 

g'S 

Number 
of  Cycles 

0.78 

23 

0.74 

27 

0.70 

50 

0.66 

120 

0.62 

280 

0.58 

600 

0.50 

2000 

0.42 

4000 

* 

0.32 

7720 

* This  value  should  be  9460  cycles  in  order  to  be  consistent 
with  the  184  ground  load  segments  which  is  16,560  cycles. 
However,  because  the  accumulative  incremental  load  factor 
falls  off  significantly  beyond  14,820  cycles  the  lower 
value  for  the  0.32  level  is  used. 


3.  Acoustic  Loads 


The  acoustic  load  levels  and  durations  are  required  to  be  established  for 
the  ground  and  flight  regimes  for  simultaneous  application  with  the  corres- 
ponding flight-by-flight  loads.  The  acoustic  spectrum  level  developed  for 
the  test  program  is  based  on  the  estimated  takeoff  spectrum  level  of  139  dB 
on  the  wing  of  a STOL  type  advanced  transport  airplane  (Reference  29).  The 
flight  spectrum  level,  however,  is  not  defined  and  is  calculated  using  the 
takeoff  noise  level  and  a variation  in  the  noise  level  with  airplane  speed 
based  on  the  fifth  power  of  the  relative  jet  exhaust  velocity. 

The  duration  of  the  acoustic  loading  is  generally  obtained  by  relating  the 
airplane  power  utilization  profile  and  the  corresponding  jet  noise  levels 
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to  an  equivalent  duration  at  the  takeoff  noise  level  (Reference  27).  For 
sonic  fatigue  analysis,  the  equivalent  duration  is  computed  using  a trade- 
off proceedure  relating  acoustic  level  to  duration  through  random  fatigue 
S/N  data  for  the  appropriate  airplane  structure.  The  power  utilization 
profile  includes  engine  maintenance  runs,  aircraft  taxi,  takeoff  and  landing 
in  addition  to  the  flight  missions  (Reference  27).  Thus,  the  progressive 
reduction  in  the  noise  level  during  the  takeoff  run  and  the  subsequent 
flight  with  increasing  airplane  speed  and  altitude  (Reference  27)  is 
accounted  for  in  the  tradeoff  procedure. 

The  duration  of  the  ground  utilization  up  to  liftoff  and  after  touchdown, 
was  estimated  at  18  seconds  per  flight  in  the  absence  of  detailed  information. 
This  represents  a total  duration  of  201  hours  for  the  40,250  flights.  The 
duration  of  the  flight  acoustic  loading  is  based  on  the  mission  profile 
data  for  a STOL  type  advanced  transport  airplane  (Reference  30).  The  flight 
noise  level  is  calculated  using  the  fifth  power  law  of  the  relative  jet  exhaust 
velocity  and  a tradeoff  proceedure  based  on  crack  propagation  data  in 
Reference  16. 

In  Reference  16  a relationship  between  the  crack  growth  rate  and  the  rms 
stress  intensity  factor  is  established  for  a panel  with  an  edge  crack 
subjected  to  random  acoustic  excitation.  Since  the  acoustic  level  is 
proportional  to  the  panel  rms  stress  and  therefore  to  the  rms  stress 
intensity  factor,  the  ratio  of  the  durations,  ti  and  t2,  of  equal  crack 
growth,  at  the  corresponding  excitation  levels,  C01  and  Gq2  (in  dB/Hz.)  is 
given  by 


lz  = 

*1 


10 


0.2  ( 


G02_  °01^ 
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Based  on  the  use  of  the  above  equation  together  with  the  flight  utilization 
data,  the  equivalent  total  flight  duration  of  130  hours  is  obtained  for  the 
40,250  flights  at  a spectrum  level  of  139  dB.  The  130  hours  represents  a 
duration  of  11.6  seconds  per  flight.  Since  the  combined  ground  and  flight 
portions  have  to  be  contained  within  28  seconds  per  flight,  the  duration 
of  the  flight  portion  is  reduced  to  10  seconds  for  the  flight  acoustic 
load  representation  in  this  program. 

The  use  of  the  accelerated  flight  spectrum  level  of  139  dB  for  simulating 
the  flight  noise  environment  essentially  assumes  that  no  interaction  effects 
on  crack  growth  are  obtained  from  the  combined  acoustic  and  flight  loads. 

The  benefits  of  a constant  flight  and  takeoff  acoustic  load  spectrum  which 
include  reduced  flight  time,  ease  of  data  interpretation,  and  simplified 
acoustic  load  simulation  outweight  the  above  limitation. 


U . Load  Tapes 


i+.l  Flight-by-Flight 

The  flight-by- flight  load  tapes  for  controlling  the  MTS  fatigue  machines 
are  generated  by  a hybrid  computer  program  based  on  the  load  histories 
developed  in  Section  III2.  The  randomization  of  the  loads  is  carried  out 
within  the  computer  program.  Each  cycle  of  the  randomization  load  is 
represented  by  a sinewave. 

The  complex  nature  of  the  combined  flight-by-flight  and  acoustic  loads 
required  the  use  of  a step-by-step  process  for  the  application  of  the 
combined  loads.  Three  flight-by- flight  load  tapes  of  increasing 
complexity  have  been  developed  for  use  in  conjunction  with  the  acoustic 
loading  for  each  set  of  the  upper  and  lower  wing  surface  loads.  The 
first  load  tape  contains  the  spectral  load  time  history,  the  second  tape 
contains  flight-by-flight  loads  but  with  the  duration  of  the  compression 
load  reduced  to  a minimum,  and  the  third  tape  contains  the  full  flight- 
by-flight  load  time  history.  The  sequence  and  level  of  each  load  peak 
is  identical  within  each  set  of  the  three  tapes.  Representative  sections 
of  the  three  load  tapes  are  illustrated  in  Figures  15  and  16  for  the  upper 
and  lower  wing  surface  loads,  respectively. 

Each  load  tape  contains  load  calibration  signal  pulses  for  counting 
flights  and  two  identical  load  tracks  for  comparison  purposes  as  a safe- 
guard against  spurious  load  pick-up  by  the  MTS  fatigue  machines . 

b.2  Acoustic 

The  acoustic  loading  is  recorded  on  one  inch  wide  magnetic  tape  for 
reproduction  by  an  Ampex  tape  recorder  that  feeds  the  signal  directly 
to  the  shaker  amplifier  independent  of  the  MTS  fatigue  machine  tape  deck. 
The  random  acoustic  loading  is  generated  by  analogue  equipment  consisting 
of  a Bruel  and  Kjaer  (B&K)  random  noise  generator,  a 3®  signal  peak  clipper 
to  prevent  shaker  overload  shut  down,  a 2 Hz  bandwidth  filter  with  a 63  Hz 
center  frequency  to  simulate  a single  degree  of  freedom  system  response, 
and  an  attenuator  to  reduce  the  signal  to  a level  acceptable  to  the  tape 
recorder.  This  method  of  generating  the  acoustic  loading  is  necessitated 
by  the  very  high  internal  damping  of  the  shaker  which  eliminates  the  panel 
resonance  peak  as  encountered  in  the  test  program.  The  resulting 
signal  history  is  typical  of  a narrow  band  random  load  history,  with  an 
appearance  of  a randomly  modulated  sinewave. 


TAPE  GA1 
SPECTRAL 


FIGURE  15.  TYPICAL  SECTION  OF  UPPER  WING  LOAD  TIME  HISTORY 


TAPE  FA1 
SPECTRAL 
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FIGURE  16.  TYPICAL  SECTION  OF  LOWER  WING  LOAD  TIME  HISTORY 


IV  CRACK  GROWTH  THEORY 


1.  Introduction 


The  analysis  and,  ultimately,  the  prediction  of  the  rate  and  path  of  crack 
growth  due  to  realistic  loading  conditions  are  of  great  importance  to  the 
study  of  failure  in  aircraft  structures.  Analysis  methods  have  been 
developed  to  predict  the  crack  growth  due  to  flight-by- flight  loads.  As 
the  problem  of  combined  acoustic  and  flight  loads  has  not  previously  been 
investigated,  no  theory  has  been  developed  which  predicts  the  effect  of 
combined  acoustic  and  flight  loads  on  crack  growth. 

In  the  absence  of  a suitable  combined  load  theory,  it  is  necessary  to 
determine  if  existing  analysis  methods  can  be  used  to  predict  the  crack 
growth  under  combined  acoustic  and  flight  loads,  and  if  not,  to  indicate 
the  direction  of  future  development.  Accordingly,  the  existing  analysis 
methods  pertaining  to  the  prediction  of  crack  growth  due  to  flight  loads 
and  acoustic  loads  are  summarized.  Equations  used  in  the  subsequent  crack 
growth  analysis  are  identified. 

2.  Flight-by-Flight  Load  Analysis 

The  programmed  crack  orientation  employed  in  all  the  combined  load  tests 
performed  in  this  study  is  in  the  most  critical  direction  for  the  flight 
loads  because  it  is  hypothesized  that  the  crack  will  be  originated  by  the 
flight-by-flight  loads.  Consequently,  the  flight-by- flight  analysis 
becomes  more  significant.  Any  significant  interaction  effects  between 
the  acoustic  and  the  flight  load  induced  panel  stresses  on  the  crack 
growth  may  be  indicated  by  deviations  from  the  theoretically  predicted 
data  beyond  normal  experimental  scatter  common  to  flight  load  crack 
growth  data. 

2.1  Crack  Growth  Laws 

Numerous  publications  on  the  subject  have  proposed  various  crack  growth 
laws  based  on  physical  models,  dimensional  considerations,  and  empirical 
observations.  References  31  and  32  offer  a detailed  review  of  the  subject. 

Based  on  the  elastic  fracture  mechanics  approach,  the  stress  intensity  factor, 
K,  is  regarded  as  the  critical  parameter  controlling  the  mechanism  of  crack 
propagation.  Of  the  various  relationships  which  define  the  rate  of  crack 
growth  as  a function  of  the  stress  intensity  factor,  the  following  crack 
growth  laws  are  the  most  popular: 
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where  da/dN  is  the  crack  growth  rate, AK  is  the  change  in  stress  intensity 
factor  of  the  load  cycle,  is  the  plane  stress  fracture  toughness, 
is  the  maximum  stress  intensity  factor  of  the  load  cycle,  (3  is  an  empirical 
constant,  and  f indicates  some  monotonically  increasing  function.  Often 
da/dN  data  is  accurately  fitted  by  a log-log  straight  line.  Accordingly, 
equation  44,  for  example,  takes  the  form 


where  Cw  and  n are  empirical  constants  and  R is  the  stress  range  ratio 
(minimum  tensile  stress  divided  by  maximum  tensile  stress). 

The  crack  growth  laws  listed  above  or  their  variations  have  been  used  with 
reasonable  success  when  employed  to  describe  crack  growth  behavior  of 
different  materials  under  constant  amplitude  load  cycles.  If,  however,  the 
load  is  of  arbitrary  varying  amplitude,  a linear  cumulative  growth  approach 

is  often  found  to  be  conservative. 

2.2  Crack  Growth  Retardation  Models 

The  phenomenon  of  crack  growth  retardation  has  been  studied  since  the  early 
1960's  and  Reference  7 offers  a significant  insight  into  the  subject. 

The  retardation  model  introduced  in  Reference  7 is  often  referred  to  as 
the  "Wheeler  Model" . 

Wheeler's  model  for  crack  growth  retardation  suggests  the  utilization  of 
a reduction  factor,  following  a high  peak  load,  on  the  subsequent  constant 
amplitude  crack  growth  rate.  The  model  can  be  described  as, 
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In  the  above  equations,  m is  an  empirical  constant  and  (re  )i  is  the  distance 
from  the  crack  tip  to  the  elastic-plastic  interface  duringFthe  jth  cycle 
after  the  application  of  a high  load.  The  parameter  rp  is  the  radius  of 
the  plastic  zone  and  is  proportional  to  the  square  of  the  ratio  of  stress 
intensity  to  tensile  yield  stress. 

= Cf-^  1+8 


(t) 


The  retardation  parameter,  Cp,  calculated  using  equation  46  accounts  for 
the  retardation  effect  N cycles  after  a high  load  is  applied. 


In  a modification  to  the  Wheeler  model,  proposed  in  Reference  33,  the  crack 
growth  rate  due  to  spectrum  loading  is  expressed  in  terms  of  the  effective 
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and  is  the  maximum  stress  intensity  factor  due  to  the  overload,  ( rp ) 

is  the  plastic  yield  zone  size  due  to  the  overload,  andSAa  is  the  sum 
of  crack  growth  since  the  overload.  In  addition,  the  coefficient  m 


(Reference  33)  by 
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where  g (AK-^h  | ,R)  is  the  function  used  in  defining  the  threshold 

stress  intensity  factor, 
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and  S is  approximately  equal  to  2.3  for  2024-T3  aluminum  and  equal  to 
2.8  for  Ti-6Al-4v. 


Another  retardation  model,  proposed  by  Willenborg,  et  al.  Reference  8, 
is  based  on  the  approach  of  an  effective  stress  intensity  factor  and  can 
be  expressed  as, 
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For  spectrum  loading  in  the  Willenborg  model,  the  crack  growth  retarda- 
tion is  expressed  in  terms  of  a stress  reduction,  a ,.  For  the  jth 
cycle  in  a load  sequence  re 


In  equation  57,  both  rep  and  rp  are  computed  for  the  current  jth  cycle. 

The  maximum  and  the  minimum  stresses  for  the  jth  cycle  (umax)j  an(^  (°min) j 
are  then  reduced  by  the  amount  (tfred)j  to  the  respective  effective  values, 
(amax)  je^  and  (fJrnin) These  effective  stress  levels  are  then  used 
together  with  appropriate  crack  growth  rate  curves  to  estimate  crack 
growth  due  to  the  jth  cycle.  By  this  means  the  Willenborg  model  becomes 
much  simplier  to  use  than  the  Wheeler  model. 

For  flights  with  a large  number  of  cycles,  where  it  is  difficult  to 
calculate  the  effect  of  retardation  cycle-by-cycle,  the  retardation  effect 
can  be  approximated  by  assuming  a constant  rer,  caused  by  the  highest  cyclic 
load.  y 


50 


This  approximation  results  in  a slight  overestimate  of  ore(j  and,  consequently, 
an  overestimation  of  the  retardation  effect.  In  other  words,  a lower  bound 
of  the  crack  growth  rate  is  obtained  using  this  assumption.  Since  the 
crack  growth  rate  calculated  with  no  retardation  provides  a corresponding 
upper  bound,  a reasonable  estimate  of  the  crack  growth  behavior  can  be 
obtained  without  an  extensive  analysis  effort. 


Using  the  constant  elastic-plastic  interface  assumption  and  selecting  the 
elastic-plastic  zone  size  produced  by  the  highest  stress,  omax,  in  the  entire 
load  spectrum  to  be  the  constant  rgp  value,  equation  53  can  be  rewritten 

as. 


j 

(0  ) = (a  ) 

red  j max  j 

1 

amax 
( °max  ^ 

_ 
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max  j I max  ^ max’ 
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where  the  function  MAX  (x,y)  takes  the  value  of  the  largest  of  its 
arguments . 


2.3  Stress  Intensity  Factor  for  Open  Hole 

The  initial  crack  in  this  study  is  programmed  into  a rivet  hole. 
Consequently,  the  stress  intensity  formula  for  an  open  hole,  given  in 
Reference  3^->  is  modified  for  an  assumed  close-fit  hole  as  discussed 
in  Reference  10.  The  resulting  Mode  I stress  intensity  factor  is  given  by: 

[tti]  63 


where  r is  the  hole  radius,  a is  the  crack  length,  and 


with  m = 1.85  and  V - —7-  > 

r+a 

or  Fq  (v)  = exp  (1.2133  -2.205 V +0.6451  V 2 ) 
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2.4  Crack  Growth  Calculation 


The  proceedure  for  calculating  the  crack  growth  due  to  flight  loads 
involves  the  computation  of  the  modified  stress  intensity  factor  (equations 
63  to  65)  corresponding  to  the  minimum  and  maximum  stress  in  each  load 
cycle.  The  corresponding  maximum  and  minimum  stress  intensity  factors 
are,  thereafter,  substituted  into  equation  44(a)  to  obtain  the  crack  growth 
rate  and  by  summation,  the  crack  growth.  With  crack  retardation  included, 
the  maximum  and  minimum  flight  load  stresses  are  replaced  in  the  above 
analyses  by  the  effective  maximum  and  minimum  stresses  given  by  equations 
59  and  60  for  the  simplified  Willenborg  retardation  model.  This  procedure 
for  calculating  the  crack  growth  has  been  computerized  employing  an 
integration  procedure  (Reference  36)  in  place  of  the  cycle-by-cycle 
summation  to  speed  up  the  computation. 

The  simplified  Willenborg  retardation  model  is  used  in  this  program  because 
of  its  ease  in  use  and  acceptable  accuracy  as  verified  later  in  this  renort. 
The  baseline  crack  growth  data,  used  in  conjunction  with  equation  44(a), 
is  obtained  from  Reference  35  for  the  7075-T6  alloy  used  in  this  program. 

The  maximum  and  minimum  flight  load  stresses  and  the  number  of  flight 
cycles  at  each  stress  level  of  the  flight  load  snectra  developed  in  section 
III2  and  used  in  the  analysis  are  summarized  in  Tables  9 and  10  for  the 
lower  and  upper  wing  surfaces,  respectively.  The  corresponding  effective 
maximum  and  minimum  stress  levels  due  to  retardation  are  included  in  the 
above  tables . 

This  analysis  indicated  a slower  than  planned  crack  growth  time  for  a 
desired  crack  growth  of  approximately  0.8  inches.  Subsequent  analysis 
indicated  that  a 60  percent  increase  in  the  stress  level  is  required  in 
order  to  complete  tests  within  a reasonable  time  period.  A comparison  of 
the  original  stress  levels  and  the  stress  levels  used  in  the  test  program, 
together  with  the  corresponding  MTS  fatigue  machine  load  levels,  are 
summarized  in  Tables  11  and  12  for  the  lower  and  upper  wing  surfaces, 
respectively.  The  test  specimen  cross  sectional  area  of  1.433  square 
inches  is  used  in  converting  the  panel  stresses  into  the  fatigue  machine 
loads . 

The  crack  growth  predictions  corresponding  to  the  increased  load  levels 
are  disucssed  and  presented  in  Section  VI,  Correlation  and  Evaluation. 

3.  Acoustic  Load  Analysis 


The  sonic  fatigue  cracks  generally  occur  along  the  rivet  line  or  the 
stiffener  edge  initiating  at  the  center  of  the  longest  panel  side.  This 
location  corresponds  to  the  peak  bending  stress  region  in  the  fundamental 
mode.  During  the  initial  crack  propagation,  the  crack  edge  is  effectively 
simply  supported, held  in  place  by  the  rough  crack  surface.  After  some 
time  when  the  rough  crack  surface  is  ground  smooth  by  the  continual 
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TABLE  9.  STRESS  LEVEIB  FOR  THE  LOWER  WING  SURFACE  SPECTRUM 
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NOTE:  The  above  spectrum  includes:  (l)  18  low  level  missions  - 72  flights 

(2)  56  typical  missions  - 112  flights 

(3)  184  G-A-G  cycles 


TABLE  10.  STRESS  LEVEIB  FOR  THE  UPPER  WING  SURFACE  SPECTRUM 
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TABLE  11.  LOWER  WING  TEST  LOAD  AND  PANEL  STRESS  DATA 


Type  of 

Original  Levels 

Final  Test  Levels 

Load 

Stress, PSI 

Load, Pound 

> 

Stress, PSI 

Load , Pound 

Mean 

Tension 

7400 

10610 

12230 

17520 

Maximum 

Varying 

5480 

7850 

9050 

12970 

Mean  + 

Maximum 

Varying 

12880 

18460 

21280 

30490 

Compression 

7500 

10750 

7500 

10750 

TABLE  12.  UPPER  WING  TEST  LOAD  AND  PANEL  STRESS  DATA 


Type  of 

Original  Levels 

Final  Test  Levels 

Load 

Stress, PSI 

Load, Pound 

Stress, PSI 

Load, Pound 

Mean 

Tension 

5000 

7170 

8000 

11460 

Maximum 

Varying 

3900 

5590 

6240 

8940 

Mean  + 

Maximum 

Varying 

8900 

12760 

14240 

20400 

Compression 

7500 

10750 

7500 

10750 
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vibration,  the  crack  edge  becomes  free  and  the  crack  will  turn  at  both  ends 
into  the  panel.  For  long  panel  edge  cracks,  relative  to  the  panel  length, 
the  crack  will  grow  into  the  panel  in  a direction  normal  to  the  panel  edge. 

With  progressive  wave  type  acoustic  excitation  parallel  to  the  long  panel 
edge  significant  stresses  can  be  obtained  in  both  the  fundamental  and  the 
(2,l)  mode  which  will  affect  the  crack  growth  rate. 

Both  the  crack  edge  condition  and  the  two  mode  panel  response  are  included 
in  the  analysis  method  presented  in  Reference  16  for  predicting  panel  edge 
crack  growth  due  to  acoustic  loading. 

3.1  Discussion  of  Analysis  Method 

The  basic  analysis  method  is  developed  for  a fixed-fixed  type  panel  with 
an  edge  crack  positioned  at  the  center  of  the  longest  panel  edge.  Com- 
posite static  beam  type  mode  shapes  are  used  to  represent  the  cracked  panel 
deflection  in  the  extended  Rayleigh  Ritz  analysis  to  determine  the  varia- 
tion of  both  the  panel  frequency  and  the  corresponding  composite  mode 
shapes  with  crack  growth. 

Expressions  relating  the  crack  tip  stress  intensity  factor  to  the  displace- 
ment at  a modal  displacement  reference  point  are  developed  from  the  com- 
posite mode  shapes.  For  the  simply  supported  crack  edge,  the  crack  opening 
mode,  Mode  I,  is  considered  to  be  dominant.  The  Mode  I stress  intensity 
factor  for  a simply  supported  crack  of  length  2a  in  the  boundary  between  a 
flexible  panel  and  a rigid  structure  was  derived  and  given  by 

K = 0.8  a x/Wa  66 

where  <x  is  the  panel  edge  cylindrical  bending  stress  remote  from  the 
the  crack  tip.  For  a panel,  the  stress,  a , is  replaced  by  the  nominal 
crack  tip  bending  stress,  aye  , as  determined  by  the  Rayleigh  Ritz  analysis. 

Random  response  theory,  similar  to  equations  17  through  20  but  in  terms  of 
displacement  instead  of  strain,  is  used  together  with  the  predicted  two 
modal  frequencies  and  the  corresponding  cracked  panel  mode  shapes  to  obtain 
the  rms  displacements  at  the  modal  reference_ points  and,  therefore,  the 
corresponding  rms  stress  intensity  factor,  K.  The  appropriate  material 
random  reverse  bending  crack  growth  data  in  form  of  the  crack  growth  rate 
equation 


da  _ F ^n  67 

dN  " L K 

which  is  used  together  with  the  rms  stress  intensity  factor  to  predict  crack 
growth  by  means  of  numerical  integrate  on. C and  n are  empirically  deter- 
mined material  constants. 

The  above  analysis  is  very  involved  ana  requires  extensive  computation. 
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Consequently,  a simplified  crack  propagation  analysis  procedure  is  suggested 
in  Reference  1 6 based  on  the  results  of  the  more  involved  analysis. 

3.2  Simplified  Analysis  Procedure 

The  method  assumes  that  the  panel  has  significant  response  only  in  the 
fundamental  mode,  which  is  a commonly  used  assumption  in  sonic  fatigue 
analysis.  The  crack  edge  is  also  assumed  to  be  simply  supported.  The 
analysis  procedure  is  as  follows: 

a.  Compute  the  maximum  panel  edge  rms  stress,  CTye , in  the  uncracked 
panel  using  the  procedure  described  in  Reference  18. 

b.  For  a crack  length  to  panel  length  ratio  B,  tbeocrack  tip  bending 
rms  stress  will  be  approximated  by  <fye(l-B2)2  and  the  correspond- 
ing rms  crack  tip  stress  intensity  factor  K given  by 

K=0.8  aye(l-B2)2\/5Fa  68 

where  a is  the  half  crack  length. 

c.  Equation  68  can  be  used  in  the  appropriate  material  crack  growth 
equation  given  by  equation  67,  resulting  in 

||  = c[0.8aye(l-B2)V^]n  69 

d.  The  crack  growth  per  stress  reversal  is  converted  to  the  crack 
growth  per  unit  time  by  the  use  of  the  fundamental  panel  frequency 
for  the  uncracked  panel.  The  cracked  panel  frequency  remains  fair- 
ly constant  for  a. value  of  B up  to  about  0.4. 

e.  The  time  taken  for.  the  crack  to  grow  from  the  flaw  size  to  any 
crack  length  up  to  40  percent  of  the  panel  length,  can  be  deter- 
mined by  numerical  integration. 

The  above  procedure  is  satisfactory  for  a value  of  B up  to  about  0.4. 

3.3  Crack  Growth  Calculation 

The  basic  random  crack  growth  data  for  reverse  bending  (equation  67)  is  not 
available  for  the  7075-T6  material.  However,  such  data  can  be 
generated  from  extensional  constant  amplitude  crack  growth  data  by  the  use 
of  the  relationship  between  extensional  and  bending  crack  growth  data 
developed  in  Reference  37;  the  Rayleigh  distribution  of  stress  peaks,  and 
linear  summation  (Miner's  Rule)  as  described  in  Reference  16.  The  result- 
ing random  crack  growth  data  over-estimates  the  Measured  crack  growth  data. 
Better  correlation  with  the  measured  data  is  obtained  in  Reference  l6  by 
the  inclusion  of  the  crack  closure  effect  (Reference  38)  in  the  analysis. 
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The  constant  amplitude  reverse  bending  crack  growth  data  was  found,  in 
Reference  l6,  to  be  highly  dependent  on  the  loading  frequency,  differing 
by  as  much  as  a factor  of  3 between  k Hz  and  175  Hz.  This  frequency  depen- 
dence requires,  by  implication,  extensional  data  of  the  appropriate  fre- 
quency to  be  developed  in  order  to  obtain  the  corresponding  random  data 
by  the  above  procedure.  Also,  the  expression  for  the  panel  edge  crack 
stress  intensity  factor  (equation  66)  may  not  be  valid  for  riveted  panels 
where  the  peak  bending  stresses  alternate  between  the  rivet  line  and  the 
frame  edge  during  each  stress  reversal. 

Due  to  the  many  sources  of  potential  error  described  above,  generation  of 
random  reverse  bending  data  using  the  procedure  in  Reference  l6  is  not 
justified.  Existing  random  crack  growth  data  can  be  used  to  provide  the 
indication  as  to  the  crack  growth  rate  for  the  7075-T6  alloy.  The  L71 
alloy  used  in  the  crack  growth  studies  in  Reference  l6  is  similar  in  com- 
position and  crack  growth  characteristics  to  the  202^-TU  alloy,  which,  in 
turn,  has  a lower  crack  growth  rate  than  the  7075-T6  alloy.  The  lower 
loading  frequency  of  63  Hz  in  this  program,  as  opposed  to  175  Hz  used  in 
Reference  l6,  should  further  increase  the  crack  growth  rate  in  the  test 
panels.  Assuming  that  the  expression  for  the  edge  crack  stress  intensity 
factor  (equation  66)  is  also  valid  for  skin  stiffener  type  structures, 
the  random  crack  growth  data  in  Reference  l6  should  represent  a lower 
bound  for  the  expected  crack  growth  from  the  acoustic  type  crack  in  Figure 
12.  The  random  crack  growth  data  indicates  that  no  crack  growth  due  to 
acoustic  loading  can  be  expected  from  the  flight  load  crack  due  to  its 
small  width,  which  can  be  considered  as  an  approximation  of  the  effective 
crack  length  for  the  acoustic  loading. 
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V TEST  PROGRAM 
1 . Introduction 


The  objectives  of  the  test  program  are  to: 

o demonstrate  the  degree  of  correlation  achieved  in  the  panel  response 
from  acoustic  and  electromagnetic  shaker  excitation, 

o correlate  the  shaker  excitation  levels  with  the  acoustic  excitation 
levels, 

o determine  the  extent  to  which  the  combination  of  acoustic  loads  with 
flight  loads  affects  crack  growth,  and 

o obtain  data  to  evaluate  current  crack  growth  analysis  methods. 

A two  phase  test  program  is  used  to  accomplish  the  above  objectives.  The 
first  phase  involves  the  demonstration  of  the  degree  of  correlation 
achieved  between  the  shaker  and  acoustic  excitation,  culminating  in  a cali- 
bration curve  between  the  acoustic  and  shaker  load  levels.  An  acoustic 
progressive  wave  tunnel,  a flexibly  mounted  electromagnetic  shaker,  and 
a single  test  panel  are  employed  in  this  phase  of  the  program. 

The  second  phase  represents  the  main  test  effort  involving  the  generation  of 
crack  growth  data  from  the  various  acoustic  and  flight  load  combinations 
to  isolate  any  effects  on  crack  growth.  Two  MTS  fatigue  machines  together 
with  two  flexibly  mounted  electromagnetic  shakers  are  used  to  apply  the 
prescribed  load  successively  to  a total  of  twelve  test  specimens. 

2.  Calibration  Tests 


The  crack  growth  rate  due  to  acoustic  loading  is  proportional  to  the  rms 
stress  (Reference  1 6)  in  the  region  of  the  crack.  The  calibration  between 
the  acoustic  and  shaker  loading  is  consequently  obtained  on  the  basis  of 
equal  rms  panel  stress  at  the  appropriate  location  on  the  test  panel.  Panel 
edge  cracks  are  used  in  the  second  phase  of  the  test  program.  It  is,  there- 
fore, necessary  to  determine  the  panel  stress  distribution  in  this  region  as 
a function  of  both  the  shaker  and  acoustic  loading.  The  rms  panel  stress  is 
also  a function  of  the  panel  frequency  and  damping  which  are,  consequently, 
of  interest. 

2.1  Instrumentation 

Six  120  ohm  strain  gages  and  one  B&K  type  8200  force  gage  are  attached  to 
the  acoustic  calibration  test  panel.  The  strain  gage  locations  are  illus- 
trated in  Figure  17.  The  B&K  type  8200  force  gage  is  attached  to  the  center 
of  the  panel.  Two  half- inch  Dynagage  type  524  microphones  are  mounted  flush 
with  the  interior  wall  of  the  progressive  wave  tunnel,  one  opposite  the 
center  of  the  test  panel  and  the  other,  which  is  used  as  a reference  microphone. 
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FIGURE  17.  CALIBRATION  PANEL  STRAIN  GAGE  LOCATIONS 
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in  an  upstream  location 


The  microphones  are  calibrated  prior  to  the  test  with  a Photocon  PC-125 
calibrator  at  lAO  dB  and  1000  Hz  frequency.  The  force  gage  is  calibrated 
by  loading  the  force  gage  with  a known  weight.  The  strain  gage  and  ampli- 
fier circuits  are  calibrated  by  unbalancing  the  bridge  with  a high  precision 
resistor  built  into  the  bridge  circuit. 

2.2  Test  Set-Up 

The  acoustic  progressive  wave  tunnel  with  a NORAIR  MK  VII  noise  generator 
and  the  electromagnetic  shaker  mounting  are  illustrated  in  Figure  18.  The 
shaker  is  suspended  by  four  tumbuckles  which  are  used  to  adjust  its  posi- 
tion relative  to  the  panel  center.  The  shaker  is  stabilized  by  two  O.O63- 
inch  thick  aluminum  flexures  attached  at  one  end  to  the  bottom  of  the 
shaker  and  at  the  other  end  to  a channel  frame  welded  to  the  progressive 
wave  tunnel.  The  flexures  and  turnbuckles  in  combination  with  the  shaker 
weight  permit  shaker  motion  only  normal  to  the  panel  surface,  therey  pre- 
venting excessive  moments  from  being  applied  to  the  force  gage.  This 
arrangement  enables  the  full  shaker  force  to  be  realized  at  the  panel  pro- 
vided that  the  shaker  drive  rod  resonant  frequency  is  considerably  greater 
than  that  of  the  panel. 

Details  of  the  acoustic  calibration  test  panel,  the  panel  fixture,  the  force 
gage  position  and  the  shaker  drive  rod  are  illustrated  in  Figure  19.  Both 
steel  and  aluminum  l/A-inch  thick  end  plates  are  used  to  bolt  the  test 
panel  and  the  frames  at  each  end  to  a stack  of  three  one-inch  thick  plywood 
boards.  The  center  section  is  cut  out  of  the  plywood  boards  to  accommodate 
the  length  of  the  frames  and  to  Just  clear  the  panel  sides.  This  method  of 
panel  attachment  provides  a good  simulation  of  the  structural  attachment 
used  later  on  in  the  program.  The  plywood  boards  are  bolted  together  and 
held  in  the  progressive  wave  tunnel  aperture  by  cantilevered  block  and 
threaded  rod  type  clamps. 

2.3  Test  Procedure,  Data  Analysis  and  Results 


The  data  needed  to  develop  a calibration  of  the  shaker  and  acoustic  excita- 
tion is  obtained  from  several  different  tests  involving  the  shaker/panel 
combination.  These  tests  involve  calibration  of  the  progressive  wave  tunnel, 
panel  damping  and  frequency,  panel  response  to  shaker  excitation, and  panel 
response  to  acoustic  excitation. 

2.3.1  Progressive  Wave  Tunnel  Calibration  - Prior  to  the  installation  of  the 
test  panel,  the  progressive  wave  tunnel  is  calibrated  with  the  test  aperture 
blocked  by  a concrete  plug  to  obtain  the  correct  tunnel  sound  pressure  levels. 
The  noise  spectrum  is  shaped  to  provide  a constant  spectrum  level, in  the 
frequency  band  of  50  Hz  to  200  Hz,  at  the  center  of  the  test  section.  The 
sound  pressure  level  is  measured  using  the  two  flush  mounted  microphones. 
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FIGURE  18.  SHAKER  AND  SPECIMEN  MOUNTING  ON  PROGRESSIVE  WAVE  TUNNEL 


FIGURE  19.  SHAKER/PANEL  DRIVING  ROD  AND  FORCE  GAGE  DETAIL 
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The  upstream  reference  microphone,  which  is  less  susceptible  to  panel 
acoustic  radiation  effects,  is  used  to  monitor  the  sound  pressure  level 
in  the  subsequent  panel  response  tests. 

The  data  is  analyzed  using  the  Hewlett-Packard  Model  8054A  third  octave 
real  time  analyzer  together  with  the  B&K  type  24l6  rms  voltmeter  to  obtain 
third  octave  sound  pressure  levels  and  the  overall  sound  pressure  level. 
Typical  spectrum  level  plots  obtained  from  the  third  octave  data  are  illus- 
trated in  Figure  20  for  both  microphone  locations. 

2.3.2  Panel  Damping  and  Frequency  - The  damping  tests  are  performed  with  the 
test  panel  installed  in  the  progressive  wave  tunnel  and  with  the  shaker 
attached  to  the  test  panel.  The  panel  is  excited  by  tapping  the  shaker  coil 
when  the  power  to  the  shaker  is  switched, in  turn,  both  on  and  off  to  produce 

a free  decay  of  the  panel  response.  The  free  decay  as  measured  by  the  strain 
gages,  is  recorded  on  magnetic  tape  and  fed  into  the  Central  Data  System 
Xerox  5 computer  for  both  a quick  look  on  the  cathode  ray  tube  (CRT)  display 
and  hard  copies  of  the  data.  Typical  free  decays  are  illustrated  in  Figure 
21  with  the  shaker  power  switched  on  and  in  Figure  22  with  the  shaker  power 
switched  off . 

The  panel  viscous  damping  coefficient  is  obtained  from  the  log-decrement  of 
the  decay.  The  measured  viscous  damping  coefficients  for  the  panel  funda- 
mental mode  are  approximately  0.4  and  0.0095  for  the  shaker  power  on  and  off 
conditions,  respectively.  The  corresponding  panel  fundamental  frequency 
obtained  by  a cycle  count  of  the  free  decay  is  approximately  70.4  Hz. 

2.3.3  Panel  Response  to  Shaker  Excitation  - A constant  amplitude  sinusoidal 
frequency  sweep  of  the  test  specimen  is  initially  performed  using  a Spectral 
Dynamics  SD  104A  sweep  oscillator  with  a SD  105A  servo  control.  The  sweep 
produced  a near  flat  panel  response  indicating  the  presence  of  high  panel 
damping. 

To  overcome  the  problem  of  high  panel  damping,  the  narrowest  available  filter 
(Spectral  Dynamics  SD  101A  tracking  filter)  with  a 2 Hz  3 dB  bandwidth  is 
used  together  with  a B&K  type  1024  random  noise  generator,  a 3cr  signal 
clipper  and  a signal  attenuator  to  generate  the  force  signal  to  drive  the 
panel.  The  filter  center  frequency  is  set  at  63  Hz. 

The  test  panel  is  excited  first  with  discrete  frequency  and  then  with  the 
narrowband  random  force,  described  above, at  typically  four  to  five  shaker 
rms  force  levels.  The  shaker  rms  force  and  the  corresponding  strain  levels 
are  measured  in  turn  with  an  rms  meter  (TSI  model  1060)  using  typically  100 
seconds  averaging  time  for  the  random  data  and  the  minimum  time  for  the 
sinusoidal  data.  A dual  beam  storage  oscilloscope  (Textronics  564)  is  used 
to  monitor  the  test  and  in  the  capture  mode,  in  conjunction  with  the  rms 
voltmeter,  to  ensure  that  the  3<J  peak  clipper  is  functioning  properly.  The 
measured  panel  edge  rms  strains  are  illustrated  in  Figure  23  as  a function 
of  the  rms  force  for  both  the  sinusoidal  and  narrow  band  random  force  exci- 
tation. 
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FIGURE  20.  TYPICAL  TUNNEL  SOUND  PRESSURE  LEVEIS  FRCM  l/3  OCTAVE  MEASUREMENTS 


FIGURE  21.  PANEL  FREE  DECAY  WITH  SHAKER  POWER  ON 


FIGURE  22.  PANEL  FREE  DECAY  WITH  SHAKER  PCWER  OFF 
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FIGURE  23.  PANEL  EDGE  STRAIN  AS  A FUNCTION  OF  SHAKER  FORCE 
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The  random  force  and  strain  data  is  recorded  on  magnetic  tape  using  a Sanga- 
mo  3568  fourteen  channel  tape  recorder.  The  force  signal  is  also  recorded 
on  magnetic  tape,  using  an  Ampex  FR  1300  tape  recorder,  for  use  in  phase  2 
tests.  Narrow  band  analysis  of  the  random  force  and  strain  data,  on  magnetic 
tape,  is  conducted  using  a Spectral  Dynamics  SD  301  Real  Time  Analyzer  in 
conjunction  with  the  SD  302  Ensemble  Averager  and  a CRT  display.  The 
analyzer  is  connected  to  the  Central  Data  System  Xerox  5 computer  for 
scaling  purposes  and  obtaining  hard  copies  of  the  data.  Typical  narrow  band 
analysis  of  the  shaker  force  (35  pounds  rms)  and  the  corresponding  panel 
edge  strain  (strain  gage  Al)  are  illustrated  in  Figures  24  and  25,  respective- 
ly* 

Peak  and  valley  probability  density  distributions  are  obtained  from  the  tape 
recorded  panel  strain  data  using  a peak  and  valley  count  routine  available 
on  the  Xerox  5 computer.  Typical  panel  edge  strain  (strain  gages  Al  and  A2) 
peak  and  valley  probability  density  distributions  with  shaker  excitation  are 
illustrated  in  Figure  2b. 

2.3.4  Panel  Response  to  Acoustic  Excitation  - The  shaker  power  is  switched 
off  for  this  test.  The  random  signal  from  the  B&K  random  noise  generator 
is  fed  through  the  General  Radio  type  1925  adjustable  level  third  octave 
filters,  to  shape  the  spectrum,  and  through  the  power  amplifier  to  the  NORAIR 
MK  VII  noise  generator.  The  sound  pressure  level  in  the  progressive  wave 
tunnel  is  measured  at  the  upstream  reference  location  using  a half- inch 
Dynagage  type  524  microphone.  The  overall  noise  level  is  measured  with  the 
B&K  type  24l6  rms  voltmeter  and  the  third  octave  levels  are  measured  with  the 
Spectral  Dynamics  real  time  analyzer.  Four  random  noise  excitation  levels, 
at  each  of  three  excitation  bandwidths  of  10,  50  and  100  Hz,  are  used  in  the 
test.  The  10  and  50  Hz  bandwidths  are  centered  at  70  Hz.  The  150  Hz  band- 
width lower  and  upper  limits  are  set  at  50  and  200  Hz,  respectively. 

The  signals  from  the  six  strain  gages  and  the  force  gage  are  fed  through 
their  respective  conditioning  systems  and  amplifiers  into  a multi  channel 
selector  switch.  The  rms  panel  strains  and  the  force  gage  outputs  are 
measured  with  the  rms  voltmeter  using  the  100  second  time  average.  The 
measured  panel  edge  rms  strains  are  presented  in  Figure  27  as  a function 
of  the  reference  microphone  overall  sound  pressure  level  for  the  three 
noise  bandwidths.  The  force  exerted  by  the  panel  on  the  shaker  coil,  as 
measured  by  the  force  gage,  is  presented  in  Figure  28  as  a function  of 
the  reference  microphone  overall  sound  pressure  level. 

The  reference  microphone,  the  force  gage  and  the  strain  gage  outputs,  which 
are  monitored  throughout  the  tests  with  a dual  beam  oscilloscope,  are  also 
recorded  onto  magnetic  tape  using  a Sangamo  tape  recorder.  Narrow  band 
analysis  of  the  magnetic  tape  data  is  performed  using  the  Spectral  Dynamics 
real  time  analyzer  and  scaled  with  the  Xerox  5 computer.  Typical  narrow 
band  analysis  of  the  reference  microphone  noise  at, 170  dB  OASPL  and  the 
corresponding  panel  strain  (strain  gage  Al)  are  illustrated  in  Figures  29 
and  30,  respectively. 
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FIGURE  2k.  TYPICAL  NARRCW  RAND  ANALYSIS  OF  SHAKER  FORCE 
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FIGURE  25.  TYPICAL  NARROW  RAND  ANALYSIS  OF  SHAKER  EXCITED  PANEL  EDGE  STRAIN 
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FIGURE  2 6.  PROBABILITY  DENSITY  DISTRIBUTION  OF  STRAIN  PEAKS  AND  VALLEYS  FOR  STRAIN  GAGES 

A1  AND  A2,  SHAKER  EXCITATION 
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FIGURE  27.  PANEL  EDGE  STRAIN  AS  A FUNCTION  OF  REFERENCE  QASPL  AND  NOISE 
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FIGURE  28.  FORCE  GAGE  OUTPUT  AS  A FUNCTION  OF  REFERENCE  QASPL  AND  NOISE 
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FIGURE  29.  NARROW  BAND  ANALYSIS  OF  ACOUSTIC  EXCITATION 
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FIGURE  30.  TYPICAL  NARRCW  BAND  ANALYSIS  OF  ACOUSTICALLY  EXCITED  PANEL  STRAIN,  STRAIN  GAGE  A1 


Typical  panel  edge  strain  peak  and  valley  probability  density  distributions 
from  acoustic  excitation  are  illustrated  in  Figure  31  for  strain  gages  A1 
and  A2.  An  estimate  of  the  panel  lateral  strain  distribution  based  on  the 
measured  strain  levels  and  strain  gage  locations  is  illustrated  in  Figure 
32.  The  strain  distribution  given  by  equation  31>  using  a reduced  panel 
width  to  reflect  the  shift  of  the  panel  edge  stress  peak  towards  the  panel 
center,  is  included  for  comparison. 

2.4  Calibration  of  Shaker  and  Acoustic  Excitation 

The  reference  microphone  overall  sound  pressure  levels,  at  selected  rms 
strain  levels  from  Figure  26,  are  first  converted  to  spectrum  levels  for 
the  three  excitation  bandwidths.  Thereafter,  the  reference  microphone 
spectrum  levels  are  corrected  to  test  section  spectrum  levels  at  the  test 
panel  resonant  frequency  using  Figure  20.  Finally,  the  rms  random  force 
levels  from  Figure  23  are  related  to  the  test  section  spectrum  level 
through  the  selected  rms  strain  levels.  The  resulting  calibration  curve  is 
illustrated  in  Figure  33* 

2. 5 Discussion  of  the  Test  Results 

2.5.1  Progressive  Wave  Tunnel  Calibration  - The  data  in  Figure  20  indicates 
that  satisfactory  spectrum  shaping  has  been  obtained  at  the  test  section 
with  the  NORAIR  MK  VII  noise  generator.  Variations  in  the  spectrum  level 

as  indicated  by  the  narrow  band  analysis  in  Figure  29  of  the  reference 
microphone  noise,  are  not  uncommon  to  acoustic  progressive  wave  test  facili- 
ties. 

The  correction  factors  obtained  from  Figure  20  appear  to  be  of  acceptable 
accuracy  considering  the  small  scatter  in  the  calculated  spectrum  levels 
for  the  three  excitation  bandwidths  illustrated  in  Figure  33* 

2.5.2  Panel  Damping  and  Frequency  - In  general,  it  is  not  possible  to  deter- 
mine the  damping  accurately,  using  the  free  decay  illustrated  in  Figure  21, 
by  the  log  decrement  method  due  to  the  limited  number  of  cycles.  The 
measured  damping  coefficient  of  0.4  can  only  be  considered  as  an  order  of 
magnitude  estimate.  Another  method  of  determining  the  damping  assumes  that 
the  frequency  of  the  decay  pulse,  in  this  instance  65  Hz,  corresponds  to  the 
frequency  of  the  much  reduced  resonance  peak.  Since  the  measured  resonant 
frequency  with  the  shaker  power  switched  off  is  70*4,  the  single-degree-of- 
freedom  theory  indicates  a damping  coefficient  of  only  0.28  for  the  5.4  Hz 
shift  in  the  damped  response  peak  frequency  (Figure  3)* 

The  measured  panel  frequency  of  70*4  Hz  without  shaker  power  is  equivalent 
to  a frequency  correction  factor  (equation  10)  of  I.33.  The  3,1  mode  was 
not  detected.  A ripple  with  a frequency  of  approximately  294  Hz  was  detected 
on  some  free  decay  traces  with  the  shaker  power  switched  on.  This  frequency 
is  considerably  higher  than  the  calculated  mode  (3,l)  frequency  and  may  be 
the  drive  rod  or  the  panel  frame  resonant  frequency. 
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FIGURE  31.  PROBABILITY  DENSITY  DISTRIBUTION  OF  STRAIN  PEAKS  AND  VALLEYS  FOR  STRAIN  GAGES 

A1  AND  A2,  ACOUSTIC  EXCITATION 
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FIGURE  32.  PANEL  STRAIN  DISTRIBUTION  RELATIVE  TO  CRACK  GROWTH  AREA 
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FIGURE  33.  CALIBRATION  CURVES  FOR  SHAKER  FORCE  AND  NOISE  EXCITATION 


The  damping  coefficient  of  0.0095,  measured  by  the  free  decay  method,  is 
in  marked  contrast  to  the  damping  coefficient  of  0.025,  obtained  from  the 
3 dB  point  bandwidth  in  Figure  30.  Resolution  errors  can  occur  if  the 
analyzer  bandwidth  is  not  sufficiently  small.  On  repeating  the  analysis  of 
the  strain  gage  response  using  an  analysis  bandwidth  of  0.3  Hz  as  opposed 
to  0.75  Hz  used  previously,  at  least  four  narrow  band  peaks  around  70  Hz 
(Figure  3k)  are  observed.  The  width  of  the  peaks  are  compatible  with  that 
expected  for  a viscous  damping  coefficient  of  about  one  percent. 

2.5.3  Shaker  and  Acoustically  Excited  Panel  Response  - The  correlation  of 
panel  edge  rms  strain  response  with  both  shaker  excitation  force  (Figure  23) 
and  acoustic  excitation  sound  pressure  level  (Figure  27)  is  generally  good 
as  characterized  by  the  small  scatter  in  the  measured  panel  edge  strains. 

The  strain  gages  A1  and  Ah,  placed  across  the  rivet  line  on  opposite  sides 
of  the  panel  (Figure  17)  in  the  location  used  as  the  reference  in  subsequent 
combined  loads  tests , produced  very  small  scatter  with  both  discrete  and 
random  shaker  excitation.  Strain  gage  Ah  failed  at  the  beginning  of  the 
acoustic  tests,  thereby  preventing  a similar  comparison  with  acoustic  ex- 
citation. Good  correlation  in  the  force  as  measured  by  the  force  gage 
during  both  the  shaker  and  acoustic  excitation  was  obtained  on  the  basis 
of  equal  rms  strain. 

The  narrow  band  analysis  (O.75  Hz  bandwidth)  of  the  panel  edge  strains 
indicates  similar  panel  response  spectra  for  both  shaker  (Figure  25)  and 
acoustic  (Figure  30)  excitation.  However,  the  strain  spectrum  level  from 
acoustic  excitation  should  be  approximately  fifty  percent  greater  than  that 
from  the  shaker  excitation  due  to  the  differences  in  the  corresponding  3 dB 
bandwidths,  since  the  rms  strain  levels  are  comparable.  The  observed  dif- 
ference in  the  spectrum  levels  is  due  to  the  presence  of  four  or  more  strain 
response  peaks  with  acoustic  excitation  (Figure  3k)  as  previously  indicated. 
The  approximate  rms  level  can  be  obtained  by  first  summing  the  peaks, 
assuming  that  the  same  bandwidth  for  each  peak  corresponds  to  a damping  co- 
efficient of  0.0095,  multiplying  by  the  bandwidth,  and  then  taking  the 
square  root.  Due  to  the  narrower  filter  and  reduced  integration  time, 
greater  statistical  errors  in  the  stress  peak  levels  are  expected.  The 
breakdown  of  a single  peak  into  multiple  peaks  can  be  caused  by  small 
geometric  differences,  for  example,  in  the  force  gage  position  or  the  rivet 
line  location  relative  to  the  frame  edge. 

The  correlation  between  the  shaker  excited  (Figure  2 6)  and  the  acoustically 
excited  (Figure  3l)  panel  edge  strain  probability  density  distributions  in- 
dicates that  a reasonable  simulation  has  been  achieved  by  the  shaker  excita- 
tion. The  distribution  of  strain  peaks  and  valleys,  illustrated  in  Figure 
31,  is  characteristics  of  riveted  panel  edge  strains  due  to  random  acoustic 
loading.  The  large  negative  strain  peaks  at  the  rivet  line  (strain  gage 
Al)  are  associated  with  the  outward  motion  of  the  panel.  The  positive 
peaks  are  limited  by  the  impingement  of  the  panel  on  the  frame  edge  during 
the  inward  motion  of  the  panel.  The  large  positive  strain  peaks  at  the 
frame  edge  (strain  gage  A2)  are  associated  with  the  bending  of  the  panel 
about  the  frame  edge  during  the  inward  motion  of  the  panel. 
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FIGURE  34.  REDUCED  BANDWIDTH  ANALYSIS  OF  ACOUSTICALLY  EXCITED  PANEL  EDGE  STRAIN 


The  estimated  lateral  panel  strain  distribution  based  on  a limited  number  of 
data  points,  illustrated  in  Figure  32,  is  influenced  at  the  panel  center  by 
the  panel  center  doublers  and  the  finite  width  of  the  force  gage  attachment. 
Consequently,  the  correlation  between  the  theoretically  assumed  and  estimated 
strain  distribution  does  not  appear  to  be  good.  However,  once  clear  of  the 
center  doubler,  the  correlation  is  expected  to  improve. 

2.5.U  Calibration  Curve  - The  calibration  curve  provides  an  acceptable  means 
for  equating  the  shaker  force  level  to  the  acoustic  excitation  level  for  use 
in  the  combined  loads  tests. 


3.  Combined  Loads  Tests 


The  specimens  used  in  the  test  program  are  designed  to  have  a constant  strain 
distribution  across  the  panel  and  to  be  free  of  bucKling  over  the  proposed 
load  range.  Verification  of  these  design  objectives  is  necessary  prior  to 
performing  the  crack  propagation  tests. 

The  crack  propagation  calculations  are  based  on  existing  constant  amplitude 
crack  propagation  data.  A baseline  crack  propagation  test  under  comparable 
loading  conditions  is  necessary  to  verify  the  use  of  the  crack  propagation 
data.  The  crack  propagation  path  and  the  variation  of  the  strain  distribu- 
tion with  crack  growth  across  the  panel  are  also  of  interest. 

The  test  specimens  are  subjected  to  variable  inplane  loads  during  the  crack 
propagation  tests.  The  panel  edge  rms  stress  produced  by  the  shaker  simu- 
lated acoustic  loading  is  affected  by  the  inplane  loads  and  has  to  be  deter- 
mined experimentally.  The  panel  fundamental  frequency  and  damping  are  re- 
quired to  correlate  the  measured  panel  edge  rms  stress  with  that  which  would 
be  obtained  from  acoustic  loading  in  order  to  indicate  the  degree  of  simu- 
lation obtained. 

Crack  propagation  tests  are  generally  performed  with  a large  number  of  test 
specimens  to  achieve  a high  degree  of  statistical  accuracy.  In  this  pre- 
liminary study,  a small  statistical  sample  represented  by  a single  test 
panel  at  each  load  condition,  is  used.  The  confidence  level,  however,  can 
be  increased,  to  some  extent,  by  the  use  of  identical  flight-by-flight  load 
cycles  and  cycle  sequences,  apart  from  the  compression  load  duration,  in 
each  set  of  three  load  tapes  for  the  upper  wing  surface  loads  and  also  for 
the  lower  wing  surface  loads.  Additional  consistency,  within  each  load  tape 
set,  can  be  obtained  by  assigning  each  load  tape  set  to  a specific  fatigue 
machine.  Variations  between  the  test  specimens  can  be  determined  by  sub- 
jecting the  test  specimens  to  constant  amplitude  cyclic  inplane  loads. 

To  isolate  the  effects  on  crack  propagation  of  the  combined  acoustic  and 
flight  loads,  the  flight-by-flight  tests  are  performed  without  and  then  with 
simulated  acoustic  loading. 
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The  use  of  two  load  levels,  represented  by  the  upper  and  lower  wing  surface 
loads,  in  conjunction  with  a constant  shaker  simulated  acoustic  loading, 
provides  additional  means,  on  a comparative  basis,  of  isolating  the  com- 
bined acoustic  and  flight  loads  effects  on  crack  growth. 

Calculations  for  crack  propagation  under  acoustic  loading  indicate  that 
crack  propagation  will  not  be  obtained  from  the  flight  load  type  crack. 
Consequently,  a typical  sonic  fatigue  crack  from  which  crack  propagation 
can  be  expected,  is  also  included  in  the  panels  which  are  subjected  only 
to  acoustic  loading, to  verify  the  forgoing  hypothesis  and  to  provide  pre- 
liminary data  on  crack  propagation  in  riveted  panels. 

A matrix  of  the  test  sequence,  test  panels,  load  tapes  and  loading  condi- 
tions employed  in  the  phase  2 test  program  is  presented  in  Table  13.  The 
increased  flight-by- flight  load  levels  used  in  the  test  program  correspond 
to  conversion  factors  of  lg  flight  load  to  12230  psl  stress  and  lg  ground 
load  to  8000  psi  stress  for  the  lower  and  upper  wing  surface  load  spectra, 
respectively.  The  compression  stress  is  maintained  at  7500  psi.  The 
calibration  signal  on  each  load  tape  is  equal  to  the  maximum  load  level. 

The  maximum  loads  are  30,^90  and  20,400  pounds  for  the  lower  and  upper  wing 
surfaces,  respectively  (Tables  11  and  12). 

The  acoustic  level  corresponding  to  the  increased  flight  loads  is  approxi- 
mately lbk  dB, which  is  equivalent  to  a shaker  force  of  35  lbs.  rms  (Figure 
33 )>  representing  a panel  rms  bending  strain  of  350  n in/in  (Figure  23)  at 
zero  inplane  load. 

3.1  Instrumentation 

Three  uniaxial  (F4,  F5  and  F6)  and  one  strain  gage  rosette  (FI,  F2  and  F3), 
illustrated  in  Figure  35>  are  employed  on  test  panel  CF3,  which  is  used  in 
the  baseline  test.  Each  strain  gage  is  wired  with  three  leads  for  tempera- 
ture compensation.  Strain  gages  F6  and  F3  are  used  to  establish  the  posi- 
tion of  the  neutral  axis.  The  strain  gage  rosette  is  used  to  define  the 
state  of  stress  in  the  test  panel  and,  consequently,  is  placed  clear  of  the 
fasteners.  Strain  gages  F3  to  F5  provide  information  on  the  strain  distri- 
bution across  the  panel.  Strain  gage  F5  is  also  used  to  indicate  the  onset 
of  buckling  should  it  occur. 

On  all  subsequent  test  panels,  strain  gage  F6  is  deleted  and  the  strain 
gage  rosette  is  relocated  with  strain  gage  FI  positioned  across  the  rivet 
line  at  the  panel  lateral  centerline,  as  illustrated  in  Figure  36.  The  FI 
strain  gage  location  corresponds  to  the  strain  gage  A1  or  Ah  location  on 
the  calibration  panel  and  is  used  to  measure  the  strain  due  to  shaker 
simulated  acoustic  loading.  The  position  of  strain  gage  F3  is  essentially 
unaltered,  whereas  strain  gage  F2  position  is  such  that  it  is  now  influenced 
by  the  fastener. 
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Table  13.  PHASE  2 TEST  PROGRAM  MATRIX 
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FIGURE  35.  BASELINE  TEST  PANEL  CF3  STRAIN  GAGE  LOCATIONS 
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FIGURE  36.  COMBINED  LOADS  AND  ACOUSTIC  LOADS  PANEL  STRAIN  GAGE  LOCATIONS 
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A single  B&K  force  gage  is  mounted  on  the  center  of  each  test  panel  to 
measure  the  shaker  simulated  acoustic  loading.  The  inplane  panel  loads 
are  measured  by  calibrated  load  cells  built  into  the  fatigue  machines. 

The  strain  gages  are  calibrated  prior  to  each  data  recording.  The  calibra- 
tion is  achieved  first  by  zeroing  the  load  cell  and  balancing  each  strain 
gage  circuit.  The  strain  gage  bridge  is  then  unbalanced  with  a high  preci- 
sion calibrated  resistor  to  provide  a known  positive  and  negative  equivalent 
strain.  Each  force  gage  is  calibrated  initially  by  loading  the  force  gage 
with  a known  weight. 

3.2  Test  Set-Up 

In  preparing  the  fatigue  machines  for  the  crack  growth  tests,  the  end  grips 
are  first  attached  to  the  fatigue  machine,  aligned  and  torqued.  The  fatigue 
machine  load  cell  is  then  zeroed  and  the  test  specimen  is  installed. 

The  installation  of  the  test  panel  in  the  fatigue  machine  involves  the  attach 
ment  of  the  panel  and  the  end  grip  attachments  to  the  end  grips  by  means  of 
the  two  1.25- inch  diameter  close  tolerance  bolts.  The  end  grip  attachments 
are  aligned  with  the  test  panel,  shimmed  where  necessary,  and  held  by  the 
5/8- inch  diameter  bolts  through  the  end  grip  attachments  and  the  panel  ends. 
The  alignment  is  facilitated  by  the  0.69“ inch  diameter  large  clearance  bolt 
holes  drilled  into  the  panel  ends.  The  test  panel  frame  web  and  lower  flange 
together  with  the  panel  skin,  are  drilled  through  the  corresponding  holes  in 
the  end  grip  attachments  and  the  fasteners  are  installed.  Four  frame  web 
stabilization  blocks,  two  to  each  end  of  the  test  specimen,  are  fitted  be- 
tween the  two  frame  webs,  shimmed  and  fastened.  The  blocks  are  located  near 
the  discontinuities  in  the  one-inch  thick  end  grip  attachment  fingers  to 
remove  web  bending  loads  introduced  by  the  discontinuities.  The  need  for 
the  stabilization  blocks  was  established  during  the  cycling  in  the  baseline 
test.  Finally,  the  1.250-inch  and  5/8-inch  diameter  bolts  are  torqued  to 
complete  the  installation.  This  complex  installation  procedure,  which  is 
dictated  by  interference  between  the  end  grip  bolts  and  the  one- inch  thick 
end  grip  attachment,  provides  good  panel  alignment.  After  the  completion 
of  the  installation,  the  load  cell  is  zeroed  in  preparation  for  strain  gage 
calibration. 

The  test  panel  installed  in  the  220,000-pound  force  MTS  fatigue  machine, 
together  with  the  control  console,  is  illustrated  in  Figure  37.  The  test 
set  up  illustrated  in  Figure  37  is  used  for  the  preliminary  compression 
tests  and  the  flight-by-flight  load  tests  which  do  not  involve  shaker  exci- 
tation. The  fatigue  machine  load  is  controlled  by  magnetic  tape  or,  for 
simpler  loads,  by  waveform  generators  in  the  control  console.  The  test  set- 
up used  in  the  combined  shaker  simulated  acoustic  and  flight-by-flight  load 
tests  is  illustrated  in  Figures  38  and  39.  The  frame  web  stabilization 
blocks  are  visible  in  Figure  38.  The  shaker  is  mounted  on  free  swinging 
tumbuckles  and  stabilized  by  two  aluminum  alloy  flexures.  The  shaker 
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FIGURE  37-  TEST  SET-TIP  FOR  BASELINE  AND  FLIGHT  LOAD  TESTS 


FIGURE  38.  TEST  SET-UP  FOR  COMBINED  ACOUSTIC  AND  FLIGHT  LOAD  TESTS 
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FIGURE  40.  TYPICAL  ARRANGEMENT  FOR  MEASURING  CRACK  GROWTH 
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load  signal  is  fed  from  the  Ampex  tape  recorder,  through  a Kronhite  hand- 
pass  filter  and  the  shaker  mobile  power  amplifier  to  the  shaker. 

The  crack  growth  is  measured  by  means  of  a travelling  microscope  and  dis- 
placement dial  arrangement,  as  illustrated  in  Figure  40,  attached  to  the 
160,000-pound  force  fatigue  machine.  The  test  specimen  and  shaker  arrange- 
ment are  also  shown  in  the  same  figure.  The  flight  loads  are  measured  by 
built  in  load  cells  and  the  shaker  force  is  measured  by  the  B&K  force  gage. 
The  load  cell,  force  gage,  and  strain  gage  signals  are  fed  through  the 
appropriate  signal  conditioning  systems  and  amplifiers  directly  into  the 
Central  Data  System  Xerox  5 computer.  The  signals  are  monitored  by  oscil- 
loscopes and  by  instant  feedback  on  the  CRT  display. 

The  load  levels  are  set  by  the  calibration  signals  on  the  flight-by-flight 
load  tape.  The  shaker  force  level  is  set  by  means  of  the  TSI  rms  meter. 

3.3  Test  Procedure 


Three  basic  test  procedures  are  used.  The  first  is  concerned  with  the  veri- 
fication of  the  panel  design  objectives  and  the  establishment  of  the  baseline 
crack  propagation  characteristics.  The  second  test  procedure  is  used  on  all 
test  panels  involving  flight-by-flight  loads  and  the  third  procedure  is  used 
with  panels  loaded  only  with  shaker  simulated  acoustic  loads.  The  test  speci- 
men response  histories  are  recorded  in  the  Central  Data  System  Xerox  5 com- 
puter. A sampling  rate  of  1500  samples  a second  is  generally  used.  For  the 
recording  of  strain,  load  cell  and  force  gage  histories,  the  sampling  rate 
is,  however,  reduced  to  750  samples  per  second.  The  flights  are  automati- 
cally counted  from  markers  on  the  load  tape. 

At  the  conclusion  of  the  test  program,  a number  of  test  specimens  are 
selected  for  fractographic  examination  of  the  crack  surface. 

3.3.1  Baseline  Tests  - The  test  procedure  used  for  the  baseline  test 
specimen  CF3  is  summarized  as  follows. 

(1)  At  zero  inplane  load,  zero  strain  gages  and  measure  precrack  length. 

(2)  Apply  compression  load  slowly  to  a panel  compression  stress  of  approxi- 
mately 7500  psi  and  back  to  zero  again. 

(3)  Monitor  strain  gage  F5  and  record  all  strain  gages  during  loading  and 
examine  quick  look  data  on  CRT  display. 

(4)  Cycle  panel  at  5 Hz  between  tension  and  compression  stress  of  89OO  and 
7500  Psi,  respectively,  and  record  strain  data. 

(5)  Examine  behavior  of  test  panel,  modify  test  panel,  as  required,  and 
repeat  above  procedure  until  a satisfactory  modification  is  obtained. 
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(6)  Cycle  panel  at  5 Hz  at  the  above  stress  levels  until  crack  growth  is 
detected,  measure  crack  length  and  note  load  cycles. 

(7)  Measure  crack  length  initially  after  every  5000  cycles  and  thereafter 
with  increasing  frequency  to  correspond  to  a crack  growth  of  approxi- 
mately 0.03  inches  between  readings  and  note  load  cycles. 

(8)  Zero  strain  gage  and  record  load  and  strain  data  at  every  half  inch  of 
crack  growth  up  to  two  inches  and,  thereafter, at  every  inch  until  the 
crack  has  propagated  across  the  panel  and  photograph  crack  progress. 

(9)  Examine  crack  surface  after  cutting  test  specimen. 

3.3.2  Flight-by-Flight  Loads  and  Combined  Loads  - The  general  test  procedure 

for  the  combined  acoustic  and  flight  loads  is  summarized  as  follows. 

(1)  Zero  inplane  load  and  all  strain  gages  and  measure  precrack  length. 

(2)  With  the  power  to  the  shaker  switched  off,  obtain  panel  free  decay  at 
every  5000  pound  inplane  load  increments  between  20,000  and  lQPOO  pounds 
of  tension  and  compression,  respectively,  by  tapping  the  shaker  coil 
and  recording  the  response  data. 

(3)  At  zero  inplane  load,  excite  test  specimen  with  shaker  using  acoustic 
load  tape  and  adjust  shaker  load  level  to  obtain  an  rms  strain  reading 
of  350  M in/in  with  a TSI  rms  meter  for  stain  gage  F3.  Note  force  level. 

(4)  Apply  inplane  load  in  increments  of  5000  pounds  from  a compression  load 
of  10,000  pounds  to  a tension  load  of  20,000  pounds,  each  time  measuring 
both  the  rms  strain  and  the  force  levels  using  a 3O  second  timg  average. 
After  test,  switch  shaker  off. 

(5)  Cycle  test  panel  at  5 Hz  between  tension  and  compression  load  levels 
of  12,760  and  10,750  pounds,  respectively,  to  initiate  crack.  Record 
the  panel  strain  response  for  all  strain  gages  for  a few  cycles. 

(6)  After  crack  initiation  measure  crack  length.  Calibrate  fatigue  machine 
with  calibration  signal  on  flight-by-flight  load  tape.  For  initial  use 
of  load  tapes,  record  load  cell  history  for  comparison  with  signal  on 
load  tape. 

(7)  Measure  crack  length  and  apply  the  shaker  and  flight-by-flight  loads. 
Record  strain  gage,  force  gage  and  load  cell  signals  and  perform  a 
peak  and  valley  count  at  the  beginning  of  the  test  for  predefined  three 
or  four  typical  flights.  Measure  crack  length  after  typically  100  to 
200  flights,  each  time  noting  the  total  number  of  flights. 

(8)  To  rewind  the  flight-by-flight  and  the  acoustic  load  tapes,  the  flight 
loads  are  switched  off  first  and  the  shaker  amplifier  gain  is  then  turned 
down. 
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(9)  When  a crack  length  of  approximately  0.8  inches  is  reached,  the  strain, 
force  gage  and  load  cell  histories  are  recorded  and  the  test  is 
stopped.  The  test  specimen  is  photographed  and  a continually  increas- 
ing static  tension  load  is  applied  to  fail  the  test  panel  and  expose 
the  crack  surface. 

When  only  flight-by-flight  load  tests  are  performed,  without  the  shaker 
being  attached  to  the  test  panel,  steps  (2),  (3),  (4)  and  (8)  are  deleted 
from  the  above  test  procedure.  The  references  to  the  shaker  load  and  the 
force  gage  are  also  to  be  deleted  from  step  (7). 

3.3.3  Shaker  Loading  - Two  test  panels  are  tested  with  shaker  simulated 
acoustic  loading.  One  of  the  test  specimens  is  simultaneously  subjected 
to  a constant  compression  load.  Since  the  two  precracks  in  the  panel  are 
at  right  angles  to  each  other  (Figures  11  and  12),  crack  initiation  by 
cycling  is  dispensed  with.  The  crack  growth  tests  with  acoustic  loading 
described  in  Reference  l6  indicates  that  rapid  crack  initiation  is  obtained 
from  a saw  cut.  The  procedure  used  is  summarized  as  follows. 

(1)  Zero  the  load  cell  and  all  strain  gages  and  measure  the  precrack 
lengths. 

(2)  Perform  free  decay  tests  at  zero  inplane  load  on  both  panels.  For 

the  panel  to  be  tested  with  a compression  load,  free  decay  tests  at  com- 
pression loads  of  5000  and  10,000  pounds  are  also  performed. 

(3)  At  zero  inplane  load,  apply  shaker  loads  from  the  acoustic  tape  to  the 
test  specimens  and  adjust  the  force  level  to  obtain  an  rms  strain  of  350 
M in/in  at  gage  F3.  For  the  test  panel  to  be  used  with  the  com- 
pression load,  also  apply  compression  loads  of  5000  and  10,000  lbs. 
while  measuring  the  rms  strain  and  force. 

(4)  At  zero  compression  lqad  for  one  test  panel  and  at  10,750  pounds 
compression  load  for  the  other  test  panel,  conduct  crack  growth  tests. 
Record  strain  gage  and  force  gage  histories  and  obtain  a peak  and 
valley  count  at  the  beginning  of  the  tests.  Note  the  crack  growth  data 
after  typically  30  to  90  minutes  of  test  time.  At  the  end  of  the  test, 
photograph  the  test  panel  and  cut  the  panel  to  examine  the  crack  surface 

3.4  Data  Analysis 

The  crack  growth  data  is  noted  during  the  test.  The  rms  strain  and  force 
levels  are  measured  with  the  TSI  rms  meter  and  also  noted  during  the  test. 

The  majority  of  the  data  is  of  the  quick  look  type  and  hard  copies  of  this 
data  are  obtained  during  or  after  the  test. 

The  Xerox  5 computer  is  used  to  analyze  the  stored  test  panel  peak  and 
valley  count  data  and  strain  response  data  for  the  flight-by-flight  and 
acoustic  load  tests.  The  analyzed  data  is  presented  as  peak  and  valley 
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distribution  and  narrow  band  strain  response  spectra.  The  narrow  band 
spectral  analysis  of  the  panel  response,  which  uses  the  Fast  Fourier  Trans- 
form, is  restricted  to  response  data  with  a constant  mean  load. 

3.5  Test  Results 


3.5.1  Baseline  Test  - Test  specimen  CF3  is  used  in  the  baseline  test.  The 
load  strain  curves  for  all  strain  gages  obtained  during  the  compression 
load  test  are  illustrated  in  Figure  4l.  Typical  load  history  and  the 
corresponding  load- strain  curve  (strain  gage  F3),  recorded  during  the  con- 
stant amplitude  cycling  at  5 Hz,  are  illustrated  in  Figures  42  and  43, 
respectively. 

The  peak  tension  and  compression  strains,  obtained  from  the  load  and  strain 
data  sampled  at  typically  inch  crack  growth  intervals  during  the  constant 
amplitude  crack  propagation  test,  are  illustrated  in  Figure  44  for  all  axial 
strain  gages  on  the  test  specimen  face  sheet.  The  ratio  of  minimum  to  maxi- 
mum peak  load  (r)  varied  between  0.854  to  0.862  during  the  constant  ampli- 
tude crack  propagation  test,  with  an  average  at  O.857. 

The  crack  growth  data  was  obtained  to  a crack  length  of  3.229  inches  after 
143,900  load  cycles.  The  test  was  stopped  after  the  test  panel  broke  from 
the  precrack  fastener  hole  to  the  panel  edge  at  a total  of  154,327  load 
cycles.  A crack  was  observed  in  the  stiffener  adjacent  to  the  panel  crack 
at  the  end  of  the  test.  The  crack  propagation  path  is  illustrated  in  Figure 
45.  The  crack  face  in  the  panel  is  initially  normal  to  the  panel  surface. 
However,  at  a distance  of  approximately  one  inch  from  the  precrack  tip,  it 
takes  on  the  appearance  of  a ridge  (or  valley)  with  the  apex  at  the  half 
panel  thickness. 

The  crack  propagation  test  data,  in  terms  of  the  crack  growth  rate  and  the 
effectiveness  stress  intensity  factor,  Keff,  is  presented  in  Figure  46. 

The  effective  stress  intensity  factor,  represented  by  the  expression  in 
the  parenthesis  in  equation  44a,  is  computed  from  the  stress  data  using 
equation  63.  Comparable  crack  propagation  data  from  Reference  35  is  in- 
cluded in  Figure  46  for  comparison. 

3.5.2  Flight-by-Flight  and  Combined  Loads  - The  flight-by-flight  and  com- 
bined load  tests  were  performed  using  panel  numbers  CF5  through  CF12.  The 
upper  wing  surface  loads  were  applied  to  panels  with  odd  numbers  using  the 
220,000  pound  force  MTS  fatigue  machine  and  the  lower  wing  surface  loads 
were  applied  to  the  panels  with  even  numbers  using  the  160,000  pound  force 
MTS  fatigue  machine. 

The  peak  tension  and  compression  stresses,  based  on  the  measured  strains 
during  the  5 Hz  crack  initiation  cycling  between  tension  and  compression 
load  levels  of  approximately  12,760  and  10,750  pounds,  respectively,  are 
illustrated  in  Figure  47.  The  crack  initiation  cycling  data  for  test  panel 
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FIGURE  1+2.  TYPICAL  LOAD  HISTORY  DURING  CYCLING 


9b 


FIGURE  43.  TYPICAL  STRAIN  LOAD  PLOT  DURING  CYCLING 


FIGURE  44.  VARIATION  OF  PEAK  AXIAL  COMPRESSION  AND  TENSILE  STRAIN  WITH 

CRACK  GROWTH  LENGTH  , PANEL  CF3 
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FIGURE  45.  CRACK  PROPAGATION  IN  TEST  SPECIMEN  CF3 

CF3  is  also  included  in  Figure  47  to  provide  additional  sampling  for 
determining  the  degree  of  variability  between  test  specimens  resulting 
from  inplane  loads. 

The  variation  of  the  panel  frequency  and  damping  with  axial  inplane  load 
obtained  from  the  free  decay  data  is  illustrated  in  Figures  48  and  49, 
respectively.  Theoretically  predicted  panel  frequency  variation  with  in- 
plane load  (equation  l4),  for  two  panel  frequency  parameter  values  of  K = 1.48 
and  1.35)  is  included  in  Figure  48.  The  variation  of  the  panel  edge  rms 
strain  with  axial  inplane  load  due  to  the  shaker  simulated  acoustic  loading 
is  presented  in  Figure  50.  All  the  panel  edge  strains  are  normalized  to  a 
panel  edge  strain  (strain  gage  Fl)  of  350  n in/in  at  zero  inplane  load 
corresponding  to  an  rms  force  level  typically  around  35  to  38  pounds. 

Typical  panel  strain  histories  due  to  the  flight  loads  without  acoustic 
excitation  (Panel  CF7)  and  with  acoustic  excitation  (Panel  CF9-2)  are  illus- 
trated in  Figures  51  and  52,  respectively,  for  strain  gage  Fl  for  the  same 
flight  on  the  load  tapes.  The  corresponding  strain  history  for  strain  gage 
F3  on  panel  CF9-2,  which  is  subjected  to  combined  acoustic  and  flight  loads, 
is  shown  in  Figure  53*  Typical  narrow  band  spectral  analyses  of  the  com- 
bined acoustic  and  flight  loads  (no  compression  loads)  strain  histories  are 
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CRACK  GROWTH  RATE.da/dN  - IN./CYCLE 


Ke„  = (I  - R)m  Kmax 

AFFDL,  FAR  73-154,  OCT  1973  R = -0.85 


STRESS  INTENSITY  FACTOR  Keff 


FIGURE  46.  BASELINE  CRACK  GROWTH  DATA,  TEST  SPECIMEN  CF3 
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FIGURE  47.  VARIATION  IN  PANEL  PEAK  STRESS  BETWEEN  TEST  SPECIMENS 
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FIGUKE  49.  PANEL  DAMPING  COEFFICIENT  AS  A FUNCTION  OF  AXIAL  INPLANE  LOAD 
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FIGUKE  51.  TYPICAL  STRAIN  HISTORY  WITHOUT  ACOUSTIC  LOADING  J STRAIN  GAGE  FI,  TEST  SPECIMEN  CF7 
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FIGURE  52.  TYPICAL  STRAIN  HISTORY  WITH  ACOUSTIC  LOADING;  STRAIN  GAGE  FI,  TEST  SPECIMEN  CF  9-2 
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FIGURE  53.  TYPICAL  STRAIN  HISTORY  WITH  ACOUSTIC  LOADING;  STRAIN  GAGE  F3,  TEST  SPECIMEN  CF9~2 


illustrated  in  Figures  5k  through  56  • 

A summary  of  the  crack  propagation  test  results  is  presented  in  Table  14. 

The  initial  crack  length  is  measured  from  the  fastener  edge  to  the  crack  tip 
after  the  crack  initiation  cycling  and  other  tests  are  completed.  Typical 
crack  propagation  patterns  are  illustrated  in  Figures  57  through  6l.  Figures 
57  (Panel  CF7)  and  58  (Panel  CFll)  represent  the  crack  patterns  for  the  upper 
wing  surface  flight  loads  without  and  with  the  acoustic  loading, respectively. 
Crack  branching  is  observed  on  panel  CFll.  The  combined  acoustic  and 
lower  wing  surface  flight  load  induced  crack  is  illustrated  in  Figure  59  for 
panel  CF12.  Figures  60  (Panel  CF5)  and  6l  (Panel  CF6)  show  the  crack  path 
for  the  combined  acoustic , tension  and  upper  and  lower  surface  spectral  loads 
(without  compression  load) , respectively.  Crack  branching  is  again  observed 
in  panel  CF5*  On  two  test  panels,  CF9-2  and  CFll,  a crack  was  initiated  at 
the  shaker  drive  hole.  The  cracks  propagated  to  the  opposite  side  of  the 
panels  to  the  precrack  as  illustrated  in  Figure  62  for  panel  CFll. 

As  indicated  in  Table  1^,  all  crack  faces  are  essentially  normal  to  the  panel 
surface  except  for  test  panels  CF8,  CF10  and  CF12  where  crack  face  incli- 
nation is  observed.  The  crack  face  inclination  in  panel  CF10,  relative  to 
the  normal  to  the  panel  surface,  started  immediately  from  the  precrack  tip. 
However,  the  crack  inclination  in  panels  CF8  and  CF12,  did  not  occur  until 
the  crack  had  propagated  to  the  stiffener  edge. 

The  crack  growth  data,  in  terms  of  the  crack  growth  rate  per  flight  as  a 
function  of  crack  length,  is  shown  in  Figure  63  for  test  panel  CF5  and 
Figure  6h  for  test  panels  CF7,  CF9-2,  and  CFll  for  the  upper  wing  surface 
loads.  The  corresponding  crack  growth  data  for  the  lower  wing  surface  loads 
is  given  in  Figure  65  for  test  panel  CF6  and  in  Figure  66  for  test  panels 
CF8,  CF10  and  CF12.  Test  panel  CF9  crack  propagation  data  is  also  included 
in  Figure  66  since  the  actual  panel  stresses  are  comparable  to  those  en- 
countered in  the  panels  subjected  to  the  lower  wing  surface  loads. 

Fatigue  cracks  occurred  also  in  the  frame  flanges  adjacent  to  the  pre- 
cracks on  all  the  test  specimens.  These  frame  cracks  are  summarized 
schematically  in  Figure  67. 

3*5.3  Simulated  Acoustic  Loads  - Two  test  specimens  were  used  in  the  simu- 
lated acoustic  crack  growth  tests  and  were  mounted  in  the  220,000  pound 
MTS  fatigue  machine.  Test  specimen  CF2  was  subjected  to  only  shaker  simu- 
lated acoustic  loading.  Test  specimen  CF^  was  subjected  to  combined  com- 
pression and  simulated  acoustic  loads. 
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NUMBER  OF  DATA  POINTS  2048 
NUMBER  OF  DATA  BLOCKS  12 
ANALYSIS  BANDWIDTH  .37  HZ 
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FIGURE  54.  NARROW  BAND  SPECTRAL  ANALYSIS  OF  STRAIN  GAGE  FI,  TEST  SPECIMEN  CF5 
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FIGURE  55.  NARROW  BAND  SPECTRAL  ANALYSIS  OF  STRAIN  GAGE  F3,  TEST  SPECIMEN  CF5 


NUMBER  OF  DATA  POINTS  2048 
NUMBER  OF  DATA  BLOCKS  12 
ANALYSIS  BANDWIDTH  .37  HZ 
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FIGURE  56.  NARROW  BAND  SPECTRAL  ANALYSIS  OF  STRAIN  GAGE  FI,  TEST  SPECIMEN  CF6 


TABLE  14.  SUMMARY  OF  FLIGHT-BY-FLIGHT  AND  COMBINED 
LOADS  CRACK  PROPAGATION  DATA 


WING 

SURFACE 

TEST 

PANEL 

INITIAL 

CRACK 

LENGTH* 

FINAL 

CRACK 

LENGTH* 

NUMBER 

OF 

FLIGHTS 

CRACK  PATH 

CF5 

0.104 

0.832 

4611 

Two  initial  cracks  with 
subsequent  branching 

UPPER 

CF7 

0.094 

0.80 

3853 

Single  Crack 

CF9-2 

0.133 

O.872 

2516 

Single  crack  on  surface, 
branching  underneath 
stiffener 

CF11 

0.126 

0.861 

2377 

Two  initial  cracks  and 
crack  branching 

CF9 

0.116 

0.805 

1145 

Single  crack 

cf6 

0.097 

0.816 

1293 

Single  crack 

LOWER 

cf8 

0.106 

O.863 

1053 

Single  crack  ** 

CF10 

0.10 

0.840 

1252 

Single  crack  ** 

CF12 

0.107 

0.912 

1006 

Single  crack  ** 

* Inches 

**  Crack  face  inclined  relative  to  panel  surface.  The  crack 
face  on  all  other  test  panels  is  approximately  normal  to 
the  panel  surface. 
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FIGURE  57.  TYPICAL  CRACK  FROM  UPPER  WING  SURFACE  FLIGHT  LOADS, TEST 

SPECIMEN  CF  7 


FIGURE  58.  TYPICAL  CRACK  FROM  COMBINED  UPPER  WING  SURFACE  FLIGHT  LOADS  AND 
SIMULATED  ACOUSTIC  LOADS, TEST  SPECIMEN  CF  11 


110 


FIGURE 


FIGURE 


59.  TYPICAL  CRACK  FROM  COMBINED  LOWER  WING  SURFACE  FLIGHT  LOADS  AND 
SIMULATED  ACOUSTIC  LOADS,  TEST  SPECIMEN  CF  12 


60.  TYPICAL  CRACK  FROM  COMBINED  UPPER  WING  SPECTRAL,  TENSION  AND 
SIMULATED  ACOUSTIC  LOADS,  TEST  SPECIMEN  CF  5 
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FIGURE  6l.  TYPICAL  CRACK  FRCM  COMBINED  LOWER  WING  SPECTRAL,  TENSION  AND 
SIMULATED  ACOUSTIC  LOADS,  TEST  SPECIMEN  CF6 


FIGURE  62.  CRACK  FROM  SHAKER  DRIVE  ROD  HOLE,  TEST  SPECIMEN  CF11 
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FIGURE  63.  MEASURED  CRACK  GROWTH  RATE  DATA,  TEST  SPECIMEN  CF5 
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FIGURE  6k . MEASURED  CRACK  GROWTH  RATE  DATA,  UPPER  WING  SURFACE  LOADS 
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FIGURE  65.  MEASURED  CRACK  GROWTH  RATE  DATA,  TEST  SPECIMEN  CF6 
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FIGURE  66.  MEASURED  CRACK  GROWTH  RATE  DATA,  LCWER  WING  SURFACE  LOADS 
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FRAME  CRACKS  ADJACENT  TO  PANEL  PRECRACKS 
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The  resonant  frequency  and  damping  coefficient  measured  from  the  free  decay 
tests  for  test  specimen  CF2  are  71.4  Hz  and  0.0069,  respectively.  Two  shaker 
force  levels,  corresponding  to  panel  rms  edge  strains  of  350  and  1+50  n in/in, 
were  used  in  this  crack  propagation  test.  The  acoustic  type  panel  edge  crack 
(Figure  68)  propagated  a total  of  0.07  inches  in  7.5  x 10°  cycles  by  the  end 
of  the  test.  The  test  was  stopped  after  a panel  center  crack  developed  from 
the  edge  of  the  center  circular  doubler  nearest  to  the  panel  side  (Figure  69). 
The  crack  propagated  approximately  parallel  to  the  panel  edge  in  both  direc- 
tions. The  panel  center  crack  propagated  much  faster  than  the  panel  edge 
crack,  both  initiating  and  reaching  a length  of  approximately  two  inches  in 
the  above  number  of  cycles.  No  crack  propagation  was  obtained  from  the  flight 
load  type  crack  normal  to  the  rivet  line.  The  rivet  line  acoustic  crack 
propagation  data  is  presented  in  Figure  70  in  terms  of  the  crack  growth  rate 
as  a function  of  the  rms  stress  intensity  factor. 

The  resonant  frequencies  of  61.5  and  40.15  Hz  and  the  corresponding  damping 
coefficients  0.007  and  0.032  were  measured  from  the  free  decays  of  panel 
CF4  at  zero  and  5000  pound  compression  load,  respectively.  The  free  decay 
at  10,750  pound  compression  load  indicated  very  high  damping.  The  panel 
edge  rms  strain  (strain  gage  Fl),  initially  set  at  350  ju  in/in  at  zero 
load,  increased  to  440  and  450  p in/in  at  compression  loads  of  5000  and  9000 
pounds,  respectively.  Unusually  large  panel  amplitudes  were  obtained  at  a 
compression  load  of  10,750  pounds.  As  a consequence,  the  compression  load 
was  reduced  until  the  large  amplitudes  abated.  Crack  growth  tests  at  a 
compression  load  of  9000  pounds  at  a panel  edge  strain  of  445  M in/in  re- 
sulted in  a slow  panel  edge  crack  growth  (Figure  JO)  from  the  acoustic  crack. 
However,  a very  rapid  panel  center  crack  growth,  in  opposite  directions 
along  the  longitudinal  centerline  of  the  panel  from  the  edge  of  the  circular 
panel  center  doublers  was  obtained.  The  panel  edge  crack  grew  0.062  inches 
while  the  panel  center  crack  both  initiated  and  grew  2.9  inches  in  0.68  x 
10°  cycles,  causing  the  panel  to  buckle.  The  crack  growth  test  was  stopped 
at  this  point.  No  crack  propagation  was  obtained  from  the  flight  load  crack. 

3.6  Discussion  of  Test  Results 

3.6.1  Baseline  Tests  - The  essentially  linear  compression  load  strain  curves 
indicated  that  panel  buckling  would  not  be  encountered  within  the  test  load 
range.  In  general,  the  load  strain  data  from  the  5 Hz  cycling  provided  more 
consistent  results  than  the  static  test. 

The  strain  data  obtained  from  the  strain  gage  rosette  (strain  gages  Fl , F2 
and  F3)  during  the  cycling  indicated  that  a near  plane  stress  state  has 
been  achieved  within  the  test  panel.  The  principal  stress  axis  is  inclined 
at  approximately  1.5  degrees  relative  to  the  longitudinal  centerline  which 
could  partially  explain  the  downward  slope  of  the  crack  propagation  path 
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FIGURE  68.  ACOUSTIC  PANEL  EDGE  CRACK,  TEST  SPECIMEN  CF2 


FIGURE  69.  PANEL  CElfPER  CRACK,  TEST  SPECIMEN  CF2 
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relative  to  the  lateral  centerline  of  approximately  four  degrees  (Figure  1+5). 
On  assuming  a plane  stress  state  and  a tension  and  compression  Young's 
modulus  of  10.3  x 106  and  10.5  x 106  psi,  respectively,  for  the  7075-T6 
material,  a test  specimen  cross  sectional  area  of  1.433  square  inches  is 
obtained  from  the  cycling  strain,  tension  load  and  compression  load  data.  The 
cross  sectional  area  from  the  panel  geometry  is  1,468  square  inches. 


Despite  the  adverse  location  of  strain  gages  F4,  just  above  the  precrack, 
and  F5,  adjacent  to  the  panel  end  doubler,  in  comparison  to  strain  gage  F3, 
only  a small  scatter  in  the  data  from  the  three  strain  gages  is  observed,  as 
illustrated  in  Figure  44  at  zero  crack  growth  length.  This  result  indicates 
that  an  even  strain  distribution  has  been  achieved  across  the  test  panel. 

With  the  exception  of  the  strain  measurements  from  strain  gage  F4,  which 
■falls  off  with  crack  growth,  the  strain  measurements  (strain  gages  F3  and  F5) 
remain  fairly  constant  up  to  a crack  length  of  approximately  two  inches. 


The  crack  growth  data  shows  good  agreement  with  the  corresponding  data  from 
Reference  35,  represented  by  the  solid  curve  in  Figure  46.  This  curve  is 
defined  by  equation  44(a)  using  the  following  values  for  the  constants  for 
the  7075-T6  sheet  material 


V 

3.44  X 10' 

m = 

0.227 

n = 

3.2 

R =- 

- 0.85 

The  crack  length  used  in  the  calculation  extends  from  the  edge  of  the 
fastener  hole  to  the  crack  tip.  A correction  factor  corresponding  to  the 
distance  between  the  fastened  hole  edge  and  the  countersunk  fastener  head 
edge  is  added  to  the  experimental  crack  length  which  is  measured  from  the 
countersunk  fastener  head  edge. 

3,6,2  Flight-by-Flight  and  Combined  Loads  - The  response  of  the  test  speci- 
men under  combined  simulated  acoustic,  inplane  and  flight  loads  is  discussed 
first. 

3. 6. 2.1  Response  of  the  Test  Specimens  — The  variation  in  the  peak  tension 
and* compression  stresses  for  all  the  test  specimens  except  for  CF9  is  quite 
small  indicating  that  good  consistency  has  been  obtained  in  the  fabrication 
and  installation  of  the  test  specimens.  The  high  stress  in  panel  CF9  is 
attributed  to  some  inadequately  tightened  fasteners  joining  the  frame  web 
to  the  fingers  on  the  end  grip  attachments,  resulting  in  elongated  frame 
web  fastener  holes.  The  elongated  frame  web  fastener  holes  caused  the 
test  specimen  face  sheet  to  be  subjected  to  a greater  portion  of  the  in- 
plane load,  thereby  increasing  the  panel  stress  by  approximately  45  percent. 
A nonlinear  hysteritic  type  load/strain  curve  is  obtained  during  the  cycling. 
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The  free  decay  response  from  panels  mounted  only  in  the  220,000-pound  force 
MTS  fatigue  machine  provided  useful  panel  frequency  and  damping  data.  The 
free  decay  data  from  the  l60, 000-pound  force  MTS  fatigue  machine  exhibited 
a beating  effect  which  prevented  accurate  analysis.  The  free  decay  data 
obtained  under  compression  load  indicates  a nonlinear  panel  response. 

The  panel  fundamental  frequency  declined  by  as  much  as  6 Hz  during  the  free 
decay  of  the  panel  response.  The  panel  frequencies,  under  compression  load 
given  in  Figure  48,  represent  averaged  frequencies  over  typically  0.3  seconds 
of  the  decay  record,  and  therefore,  overestimate  the  small  amplitude  panel 
frequency  by  a small  amount.  No  nonlinear  effects  are  observed  when  the  panel 
is  under  a tension  load. 

The  zero  inplane  load  panel  frequencies  are  similar  to  that  measured  on  the 
calibration  panel,  and  confirm  that  the  empirical  frequency  correction 
factor  K is  approximately  1.33*  A comuarison  of  the  measured  and  theoreti- 
cally predicted  variation  (K  = 1.35)  in  the  edge  stiffened  panel  fundamental 
frequency  indicates  that  a reasonable  estimate  of  the  panel  fundamental  fre- 
quency can  be  obtained  using  equation  14  with  the  correct  frequency  correc- 
tion factor. 

The  relatively  low  panel  damping  coefficient  of  O.OO63,  for  zero  and  tension 
inplane  loads  (Figure  49),  in  comparison  to  the  damping  coefficient  of  0.0095, 
measured  for  the  calibration  panel  in  the  progressive  wave  tunnel,  may  be  due 
to  differences  in  acoustic  radiation  effects  between  the  two  test  set-ups. 
Similar  trends  were  obtained  in  damping  tests  on  honeycomb  panels  in  Refer- 
ence 39,  where  considerable  reductions  in  panel  damping  was  obtained  by  the 
removal  of  large  baffles  (simulating  large  cross  section  acoustic  progressive 
wave  tunnels)  from  around  the  test  panels.  The  small  change  in  the  panel 
damping  coefficient  will  not,  however,  affect  the  shaker  excited  panel  re- 
sponse in  view  of  the  simulation  used. 

The  use  of  strain  gage  FI  rms  level  for  establishing  the  shaker  force  levels 
for  subsequent  simulated  acoustic,  inplane  load,  and  flight-by-flight  load 
tests  is  selected  to  ensure  a consistent  rms  strain  level  around  350  n in/in 
for  all  the  test  specimens  at  zero  inplane  loads.  This  procedure  is  further 
justified  by  the  small  scatter  between  the  panel  edge  gages  A1  and  a4  on  the 
calibration  test  panel.  The  resulting  panel  rms  force  level  is,  on  the 
average,  higher  than  that  measured  on  the  calibration  panel.  The  three  force 
gages  failed  during  the  testing  from  unknown  causes.  Free  decay  and  shaker 
excited  strain  peaks  indicated  force  levels  well  within  the  capability  of 
the  force  gage. 

The  variation  in  the  panel  rms  strain  with  inplane  load,  as  illustrated  in 
Figure  50,  is  similar  to  that  predicted  by  linear  theory  for  acoustic  exci- 
tation. Equation  19  is  used  together  with  the  averaged  frequency  data 
(Figure  48)  and  the  damping  data  (Figure  49),  to  which  an  increment  of  0.0032 
has  been  added  to  represent  the  panel  damping  in  the  progressive  wave  tunnel, 
to  predict  the  variation  in  panel  rms  stress  with  inplane  load  as  illustrated 
in  Figure  50*  The  predicted  variation  is  normalized  at  an  rms  strain  of  350 
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H in/in  at  zero  inplane  load.  It  is  observed  from  Figure  50  that  during  the 
application  of  the  shaker  simulated  acoustic  and  the  ground  and  flight  spec- 
tral loads,  the  panel  edge  rms  strain  from  the  shaker  excitation  varies 
approximately  between  240  and  380  /i  in/in  and  between  200  and  380  m in/in  for 
the  upper  and  lower  wing  surface  spectral  loads,  respectively. 

The  typical  axial  strain  history  (strain  gage  F3)  under  combined  acoustic 
and  flight  loading  indicates  that  the  panel  edge  axial  strain  is  not 
affected  by  acoustic  loading  as  predicted  by  equation  15 . Comparison  of 
the  lateral  strain  history  (strain  gage  Fl)  without  and  with  the  simulated 
acoustic  loading  in  Figures  51  and  52,  respectively,  indicates  the  effect 
of  the  acoustic  loading  on  the  lateral  panel  edge  strain.  A similar  compari- 
son can  be  made  from  the  corresponding  narrow  band  strain  response  spectra 
illustrated  in  Figures  54  through  5 6*  Although  nearly  25,000  samples  are 
used  in  the  analysis,  the  samples  represent  only  three  flights  and  thus 
provide  only  a limited  representation  of  the  flight  loads.  The  strain  data 
from  panels  CF5  and  CF6  are  used  in  the  spectral  analysis  to  avoid  changes 
in  the  mean  strain  level  present  in  the  response  of  the  remaining  test  panels. 

3. 6. 2. 2 . Upper  Wing  Surface  Flight  Load  Produced  Crack  Details  - The  crack  propa- 
gation path  for  panel  CF7,  which  is  subjected  only  to  flight  loads  is  very 
straight  as  illustrated  in  Figure  57*  An  inplane  tension  load  was  applied  to 
the  panel  to  open  the  crack  for  the  photograph.  The  crack  face  is  normal  to 
the  panel  surface.  No  fatigue  markings  are  visible  on  the  photograph  of  the 
crack  face  (Figure  71)  with  a 5 times  magnification.  Examination  of  the 
crack  surface  by  means  of  an  electron  microscope  with  a magnification  of 
65OO  times  reveals  distinct  markings  corresponding  to  the  mean  load  changes 
between  each  flight  as  illustrated  in  the  fractographs  in  Figure  71.  The 
fractographs  are  taken  from  two  locations  halfway  between  the  precrack  tip 
and  the  frame  edge.  The  crack  growth  rate  of  approximately  80  x 10“ 6 inches 
per  flight,  measured  from  the  fractographs,  corresponds  well  with  the  crack 
growth  rate  of  100  x 10"°  inches  per  flight  obtained  from  Figure  64  at  a 
crack  length  of  0.25  inches. 

The  surface  crack  propagation  path  for  test  panel  CF9-2,  which  is  subjected 
to  simultaneous  flight  and  simulated  acoustic  loading  during  the  take-off 
segment  of  each  flight,  is  straight  but  inclined  slightly  downward  similar 
to  the  initial  crack  growth  path  on  panel  CF3.  Examination  of  the  crack 
surface  revealed  that  the  crack  branched  on  the  panel  surface  below  the 
stiffener.  The  crack  branching  location  is  visible  in  the  magnified  photo- 
graph of  the  crack  surface  in  Figure  72.  The  branching  crack  extends  down 
to  almost  the  half  thickness  of  the  face  sheet.  The  fractographs  of  the 
crack  surface  indicate  similar  markings  to  those  observed  on  panel  CF7* 

The  crack  which  developed  on  the  opposite  panel  side  to  the  precrack  at 
the  panel  center  drive  rod  bolt  hole  was  stop  drilled  at  a length  of  one 
inch  and  twice  after  that  but  continued  to  grow.  It  did  not  reach  the 
strain  gage  rosette  by  the  end  of  the  test.  In  view  of  the  relatively  con- 
stant strain  distribution  up  to  a crack  length  of  2 inches,  measured  on 
panel  CF3  (Figure  44),  the  panel  center  crack  most  likely  did  not  influence 
the  crack  propagation  from  the  precrack. 
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TYPICAL  FRACTOGRAPHS  OF  CRACK  SURFACE  MAGNIFIED  65OOX 
FIGURE  71.  CRACK  SURFACE  AND  FRACTOGRAPHIC  DETAIL,  TEST  SPECIMEN  CF7 
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TYPICAL  FRACTOGRAPHS  OF  CRACK  SURFACE  MAGNIFIED  65OOX 


FIGURE  72.  CRACK  SURFACE  AND  FRACTOGRAPHIC  DETAIL,  TEST  SPECIMEN  CF9"2 


125 


Two  cracks  initiated  from  the  corners  of  the  precrack  in  panel  CF11.  One  of 
these  cracks  continued  to  propagate  at  a large  angle  to  the  panel  lateral 
centerline  until  branching  occurred  (Figure  58)  after  which  a straight  crack 
propagation,  with  a small  incline,  developed.  This  test  panel  was  subjected 
to  simultaneous  acoustic  and  flight  loading  encompassing  both  the  flight 
and  ground  segments  in  each  flight.  A second  crack  (Figure  62),  similar 
to  the  one  in  panel  CF  9-2,  developed  from  the  panel  center  shaker  drive 
rod  bolt  hole.  This  crack,  with  a length  of  approximately  1.4  inches,  most 
likely  did  not  influence  the  precrack  propagation  for  the  reasons  previous- 
ly discussed.  After  the  discovery  of  the  second  crack,  the  circular  doublers 
were  removed  from  the  remaining  test  panels,  the  center  hole  was  reamed  and 
deburred,  and  the  doublers  were  rebonded  after  cleaning  the  surfaces. 

Test  panel  CF5  was  subjected  to  the  same  simulated  acoustic  and  spectral 
loads  as  panels  CF9-2  and  CF-11  except  for  the  compression  load,  which  was 
deleted.  Double  crack  branching  was  obtained  (Figure  60)  before  the  crack 
assumed  a reasonably  straight  but  slightly  inclined  crack  growth.  The  two 
crack  branch  locations  are  visible  in  the  magnified  crack  face  photograph  in 
Figure  73.  The  markings  on  the  fractographs  are  generally  uniform  over  the 
surface  and  correspond  to  the  high  flight  load  spectral  peaks. 

Test  panel  CF9  encountered  flight  load  induced  stresses  which  are  approximate- 
ly 45  percent  greater  than  the  flight  load  stresses  in  the  test  panels  pre- 
viously discussed  and,  therefore,  are  at  a comparable  level  with  the  lower 

wing  surface  flight  load  induced  stresses.  A straight  but  slightly  inclined 
crack  growtn  path  was  ootained.  The  crack  face  is  only  slightly  inclined 
to  the  normal  to  the  panel  surface. 

3. 6. 2. 3 Lower  Wing  Surface  Flight  Load  Produced  Crack  Details  - The  crack  pro- 
papgation  path  for  panel  CF8  without  simulated  acoustic  loading  is  similar  to 
panel  CF12  (Figure  59)  which  is  subjected  to  simulated  acoustic  loading  for 
both  the  flight  and  ground  segment  in  each  flight.  The  crack  face  in  both 
panels  is  initially  normal  to  the  panel  surface  but  rotated  to  an  angle  of 
45  degrees  near  the  frame  edge.  The  point  of  rotation  is  indicated  in 
Figure  59  hy  the  change  in  crack  growth  direction.  It  is  also  visible  in 
the  magnified  photograph  of  the  crack  face  in  Figure  74.  The  fractographs 
in  Figure  74  also  show  crack  growth  marks.  The  markings  which  differentiate 
the  individual  flights  are  difficult  to  identify  due  to  the  high  crack 
growth  rate  associated  with  the  panel. 

The  crack  path  for  panel  CF10  is  a shallow  curve  associated  with  the  initia- 
tion of  the  crack  face  rotation  to  a 45  degree  angle  relative  to  the  panel 
surface  directly  from  the  precrack.  The  crack  direction  is  essentially 
normal  to  the  rivet  line. 

The  crack  path  for  panel  CF6,  illustrated  in  Figure  6l,  is  essentially 
straight  with  a slight  incline  in  the  same  direction  as  all  the  test  panels 
including  panel  CF3.  The  crack  face  is  normal  to  the  panel  surface. 


126 


V\V. 


FRAME  CRACK 


FASTENER 


HOLE 


END  OF  CRACK 

■growth* 


CRACK  SURFACE  MAGNIFIED  5X 


1/i 

TYPICAL  FRACTOGRAPHS  OF  CRACK  SURFACE  MAGNIFIED  65OOX 
FIGURE  73.  CRACK  SURFACE  AND  FRACTOGRAPHIC  DETAIL,  TEST  SPECIMEN  CF5 
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FIGURE  T1*-  CRACK  SURFACE  AND  FRACTOGRAPHIC  DETAIL,  TEST  SPECIMEN  CF12 


128 


3.6.2.U  Measured  Crack  Propagation  Rates  - The  crack  propagation  rate  data  is 
calculated  from  measured  crack  propagation  data  at  typically  every  third 
measured  data  point  to  smooth  the  crack  growth  curves  and  is  presented  as  a 
function  of  the  crack  length  from  the  fastener  hole.  A least  squares  curve 
fit  is  used  to  establish  a curve  through  the  calculated  crack  growth  rate  data 
points . 

The  crack  growth  rate  data  for  the  upper  wing  surface  test  panels  CF7,  CF9-2 
and  CF11,  combined  in  Figure  6U,  indicate  a slight  increase  in  crack  growth 
rate  for  panels  CF9-2  and  CF11,  which  were  subjected  to  simultaneous  flight 
and  simulated  acoustic  loading.  The  increase  in  crack  growth  rate,  however, 
is  very  small  and  within  the  expected  data  scatter  band.  Consequently,  no 
significant  effect  on  the  crack  growth  rate  can  be  attributed  to  the  combined 
acoustic  and  flight  loads.  The  crack  growth  rate  data  from  panel  CF5  (Fig- 
ure 63),  without  the  compression  load  cycle,  is  also  comparable  to  the  above 
panels  up  to  a crack  length  of  0.6  inches,  after  which  the  crack  growth  rate 
increases  at  a slower  rate. 

The  crack  growth  rate  for  the  upper  wing  surface  flight  loads  for  test  panels 
CF8,  CF10  and  CF12  together  with  test  panel  CF9  are  summarized  in  Figure  66. 

Two  distinct  crack  propagation  rates  are  obtained.  The  two  test  panels  CF8 
and  CF12,  in  which  the  crack  rotated  to  a U5  degree  angle  on  passing  the  frame 
edge  are  comparable  and  exhibit  an  accelerated  crack  growth  rate  relative  to  the 
cracks  in  panels  CF9  and  CF10.  The  crack  growth  rates  for  test  panels  CF9,  CF10 
and  CF6,  Figure  6 5,  are  also  comparable.  The  crack  growth  data  indicates  that 
the  combined  simulated  acoustic  and  flight  loads  have  no  effect  on  the  crack 
growth  in  these  test  panels. 

3.6.2. 5 Crack  Branching  - Crack  branching  is  present  only  in  the  panels  which 
are  subjected  to  the  combined  acoustic  and  upper  wing  surface  flight  loads. 

No  crack  branching  is  observed  on  panel  CF7  which  is  subjected  only  to  the 
upper  wing  surface  flight  loads.  Consequently,  it  is  deduced  that  the  crack 
branching  is  caused  by  the  acoustic  loading  in  combination  with  the  flight 
loading.  Additional  supporting  evidence  is  obtained  by  considering  the  varia- 
tion of  the  stress  field,  produced  by  the  combined  loading,  through  the  thick- 
ness of  the  panel.  The  acoustic  loading  produces  a stress  field  at  the  panel 
edge,  in  the  region  of  the  precrack,  which  is  a maximum  at  the  surfaces  of  the 
panel  (either  tension  or  compression)  and  zero  near  the  half  thickness  of  the 
panel.  Thus,  the  combined  acoustic  and  flight  loads  produce  a biaxial  stress 
field  at  the  surface  of  the  panel  (equations  29  and  30)  and  a uniaxial  stress 
field  at  the  half  thickness  of  the  panel.  Any  acoustic  load  effects  on  the 
crack  growth  are,  therefore,  expected  to  decrease  progressively  through  the 
thickness  of  the  panel  and  stop  near  the  panel  half  thickness.  The  crack 
branching  penetration  depth  in  Figures  72  and  73  appears  to  corroborate  the 
above  reasoning.  Crack  branching  also  produces  a greater  scatter  in  the  crack 
growth  data  as  evident  in  Figure  63. 

No  crack  branching  is  observed  on  the  panels  subjected  to  the  combined  acous- 
tic and  the  much  greater  (see  Tables  11  and  12)  lower  wing  surface  tension 
flight  loads.  Consequently,  the  ratio  of  the  tension  flight  stress  levels 
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to  the  acoustic  stress  levels  is  a significant  parameter  in  determining  the 
onset  of  crack  branching. 

3. 6. 2. 6 Frame  Fretting  Crack  Details  - Fretting  initiated  fatigue  cracks  in 
the  panel  frames  are  drawn  in  Figure  67,  magnified  approximately  five  times. 
These  fretting  cracks  were  all  initiated  at  the  fastener  hole  in  the  frame 
adjacent  to  the  precrack  root.  A typical  fretting  initiated  crack  is  visible 
in  the  magnified  crack  face  photograph  in  Figure  73.  No  frame  crack  initia- 
tion and  crack  growth  data  could  be  obtained  since  the  crack  occurred  in  an 
inaccessible  location. 

3.6.3  Simulated  Acoustic  Loads  - The  frequency  and  damping  coefficient  mea- 
sured on  test  specimen  CF2  at  zero  inplane  loads  are  of  comparable  level  to 
those  measured  on  the  other  test  specimens.  The  resonant  frequency  of  61.5  Hz 
at  an  indicated  zero  inplane  load  measured  on  panel  CF4  is,  however,  too  low. 
The  behavior  of  the  test  specimen  CF4  under  compression  load,  both  during  the 
free  decay  and  the  shaker  excitation  tests,  indicated  the  presence  of  panel 
buckling. 

Examination  of  the  test  specimen  CF4  free  decay  response  at  supposedly  zero 
inplane  load,  after  the  completion  of  the  test  program,  revealed  the  presence 
of  a nonlinear  panel  response,  associated  with  compression  load,  indicating 
that  incorrect  load  readings  were  obtained  from  the  load  cell  circuit.  The 
frequency  of  6l.5  Hz  and  a damping  coefficient  of  0.007  are  consistent  with 
a compression  load  of  approximately  10,000  pounds  based  on  the  test  data  in 
Figures  U8  and  1+9. 

The  crack  growth  data  from  test  specimen  CF2,  at  a panel  edge  rms  strain 
level  of  350  p in/in  from  only  the  shaker  simulated  acoustic  loading  (Fig- 
ure 70),  indicated  a rapid  fall  off  in  the  crack  growth  rate  until  the  crack 
appeared  almost  to  have  stopped  growing  after  a total  of  6 x 10^  random  load 
cycles.  This  period  of  crack  growth  is  associated  with  the  change  in  the 
crack  front  from  the  vertical  precrack  to  a bending  type  crack  involving 
high  surface  crack  growth  and  low  growth  at  panel  half  thickness.  Increas- 
ing the  panel  edge  rms  strain  level  to  1+50  \x  in/in  produced  the  desired 
increase  in  cracic  growth  rate.  The  test  was  stopped  after  a total  of  7.7  x 
106  random  load  cycles  when  a vertical  crack  was  discovered  at  the  edge  of 
the  panel  center  circular  doubler.  The  resulting  crack  growth  rate 
in  terms  of  the  panel  edge  stress  intensity  factor  (equation  66)  is  a 
factor  of  three  less  than  that  obtained  from  a panel  edge  crack  along  a 
fixed  panel  edge  (Reference  l6). 

The  slow  crack  growth  relative  to  that  reported  in  Reference  l6  is  more 
suprising  as  the  L71  alloy  used  in  that  study  is  comparable  to  2024-T4 
alloy,  which  has  a slower  crack  growth  rate  than  the  7075-T6  material.  The 
fatigue  life  of  the  shaker  attachment  was  estimated  at  2 x 10T  cycles  at  a 
stress  level  of  4500  psi  based  on  the  random  SN  data  in  Reference  19.  The 
unexpected  slow  crack  growth  rate  enabled  the  load  cycles  to  approach  the 
design  life  of  the  shaker  attachment. 
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The  crack  front  in  test  specimen  CF4,  at  an  indicated  compression  load  level 
of  9000  pounds  and  a rms  strain  level  of  41+5  M in/in,  did  not  complete  the 
transition  to  a tending  crack  front  (Figure  TO)  when  a vertical  crack  developed 
from  one  lateral  edge  of  the  round  panel  center  doublers  after  a total  of 
0.7  x 10°  load  cycles.  The  crack  grew  to  a length  of  2.9  inches  in  the  last 
half  hour  of  testing,  causing  the  panel  to  buckle. 
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VI.  CORRELATION  AND  EVALUATION 


1.  Introduction 


The  complex  problem  of  theoretically  predicting  the  stiffened  panel  re- 
sponse to  shaker  simulated  acoustic  excitation  and  simultaneously  applied 
axial  inplane  loads  resulted  in  the  use  of  a semi-empirical  method,  based 
on  empirical  data  from  stiffened  panels,  in  the  design  of  the  test  speci- 
mens. The  accuracy  of  the  semi-empirical  method  is  evaluated  against  test 
data. 

In  addition,  this  section  covers  the  most  significant  objective  of  compar- 
ing the  existing  crack  growth  analysis  methods  with  the  measured  crack  growth 
data.  Flight-by-flight  load  crack  growth  analysis,  as  described  in  Section 
IV  2.4,  has  been  performed  on  all  panels  subjected  to  flight  loads  and  combined 
loads,  both  without  cracK  retardation  and  using  the  Willenborg  crack  retard- 
ation model.  The  crack  growth  due  to  simulated  acoustic  loading  is 
evaluated  on  the  basis  of  available  test  data  for  other  aluminum  alloys  be- 
cause of  the  lack  of  baseline  crack  growth  data  at  comparable  loading  fre- 
quency for  the  7075-T 6 material. 


2.  Panel  Response 


As  indicated  in  the  discussion  of  the  test  results  in  the  previous  section, 
equation  14  provides  a reasonable  estimate  of  the  variation  of  panel  funda- 
mental frequency  with  inplane  load  (Figure  48)  when  using  the  correct  value 
for  the  panel  edge  frequency  parameter,  K.  The  relatively  low  value  of  1*33 
for  the  frequency  parameter,  as  opposed  to  the  more  typical  value  for  stiffen- 
ed airplane  type  panels  of  1.48,  indicates  an  increased  panel  edge  flexibil- 
ity probably  associated  with  frame  web  bending. 

The  electromagnetic  shaker  simulation  used  in  this  program  did  not  repro- 
duce the  variation  of  panel  resonant  frequency  with  inplane  load  that 
would  have  occurred  with  acoustic  excitation.  This  variation  in  the  panel 
frequency,  apart  from  any  modal  curvature  changes  which  would  affect  the 
stress  distribution  in  the  panel,  is  responsible  together  with  the  panel 
damping,  for  the  variation  in  the  panel  rms  stress  with  inplane  load.  In 
spite  of  this  limitation,  good  correlation  in  the  shape  of  the  measured 
panel  edge  rms  stress  variation  with  inplane  load,  when  compared  with  the 
theoretically  predicted  variation  (Figure  50),  has  been  obtained.  Thus, 
the  most  critical  parameter  affecting  crack  propagation,  namely  the  stress 
level,  has  been  adequately  simulated  over  the  range  of  inplane  load  used 
in  this  program. 

The  correlation  between  the  measured  panel  edge  stress  and  the  predicted 
panel  edge  stress,  using  equations  27  and  the  skin-stringer  design  nomo- 
graph in  Reference  19,  due  to  acoustic  loading,  is  complicated  by  the 
presence  of  at  least  four  peaks  in  the  panel  response  (Figure  34)  around 


132 


69  Hz.  On  assuming  a single  mode  response  normally  obtained,  equation  27 
in  conjunction  with  the  design  nomographs  over-estimate  the  rms  panel  edge 
stress  by  4l  percent  for  a damping  coefficient  of  0.0095.  This  difference 
of  4l  percent  is  within  the  accuracy  of  current  sonic  fatigue  analysis 
methods . 

The  correlation  curve  between  the  shaker  force  and  the  acoustic  spectrum 
level  in  Figure  33  is  obtained  on  the  basis  of  equal  panel  edge  rms  strain. 
The  mean  square  panel  edge  stress  due  to  the  shaker  force  and  also  the 
acoustic  force  are  given  theoretically  by  equation  25  when  the  generalized 
force  term,  GFj^,  is  replaced  by  the  corresponding  shaker  and  acoustic 
generalized  force  terms  given  by  equations  37  and  35,  respectively,  and 
the  corresponding  frequencies  and  damping  coefficients  are  used.  Since  a 
2 Hz  bandwidth  shaker  excitation  is  used,  the  expression  for  the  mean 
square  panel  edge  stress  (equation  25),  due  to  shaker  excitation,  must  be 
multiplied  by  a bandwidth  correction  factor  obtained  from  Figure  4.  These 
expressions  for  the  mean  square  panel  edge  stress  can  be  rearranged  on 
the  basis  of  equal  mean  square  panel  edge  stress  to  produce  the  ratio  of 
shaker  force  level  to  acoustic  spectrum  level  for  comparison  with  the  cali- 
bration curve  in  Figure  33*  With  shaker  power  on  and  power  off  damping 
coefficients  of  0.4  and  0.0095,  respectively,  the  above  procedure 
over-estimates  the  ratio  of  the  electromagnetic  shaker  force 
level  to  the  acoustic  spectrum  level  by  approximately  37  percent. 

The  degree  of  correlation  is  reasonable  in  view  of  the  uncertainty  as  to 
the  actual  value  of  the  damping  coefficient  measured  with  the  shaker  power 
switched  on. 


3.  Crack  Growth  Prediction 


The  expression  for  the  stress  intensity  factor,  given  by  equation  63,  is 
used  together  with  the  baseline  crack  growth  data  to  obtain  the  crack 
growth  rate  in  terms  of  the  effective  stress  intensity  factor.  Good  corre- 
lation has  been  obtained  between  the  baseline  crack  growth  rate  data  and 
comparable  published  data  (Reference  35),  as  illustrated  in  Figure  46. 

For  the  flight-by-flight  loads  and  the  combined  loads,  the  results  of  the 
crack  growth  analysis  procedure  discussed  in  Section  TV  2.4  are  presented 
graphically  in  terms  of  the  total  number  of  flights  ana  the  total  crack 
length.  The  results  of  the  analyses  are  illustrated  in  Figures  75  and  76 
for  the  upper  wing  surface  loads  and  in  Figures  77  and  78  for  the  lower 
wing  surface  loads  together  with  the  corresponding  crack  growth  data. 

The  agreement  between  the  measured  and  predicted  crack  growth  without  re- 
tardation for  test  specimen  CF5,  subjected  to  simulated  acoustic  loads, 
and  the  upper  wing  surface  spectral  loads,  without  the  compression  loads, 
is  reasonable  since  crack  retardation  effects  are  generally  not  signifi- 
cant for  this  type  spectral  loading.  In  general,  good  agreement  between 
the  measured  and  predicted  crack  growth  has  been  obtained  for  the  test  speci- 
mens CF5,  CF7 , CF9-2  and  CF11  which  are  subjected  to  the  upper  wing  surface 
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FIGURE  75.  COMPARISON  OF  MEASURED  AND  CALCULATED  CRACK  GROWTH,  TEST  SPECIMEN  CF5 


z 

I 

X 
H 
< J 
z 


CJ 

< 

cc 

CJ 


135 


FIGURE  76.  COMPARISON  OF  MEASURED  AND  CALCULATED  CRACK  GROWTH,  UPPER  WING  SURFACE  LOADS 
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FIGURE  77.  COMPARISON  OF  MEASURED  AND  CALCULATED  CRACK  GROWTH,  TEST  SIECIMEN  CF6 
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FIGURE  78.  COMPARISON  OF  MEASURED  AND  CALCULATED  CRACK  GROWTH,  LOWER  WING  SURFACE  LOADS 


flight-by-flight  loads.  The  measured  crack  growth  data  for  test  specimen  CF7, 
which  is  subjected  only  to  flight  loads,  falls  between  the  crack  growth  pre- 
dicted by  the  Willenborg  retardation  model  and  the  crack  growth  predicted  with- 
out crack  retardation.  For  test  panels  CF5,  CF9-2  and  CF11,  which  are  subjec- 
ted to  the  combined  acoustic  and  upper  wing  surface  loads,  good  agreement 
between  the  measured  and  predicted  crack  growth  without  crack  retardation  is 
obtained.  Only  the  axial  flight  load  induced  stresses  are  included  in  these 
calculations.  The  generally  good  agreement  between  the  measured  and  predicted 
crack  growth  data  for  all  the  test  specimens  subjected  to  the  upper  wing  sur- 
face flight  loads,  indicates  that  no  significant  effects  on  crack  growth  can 
be  attributed  to  the  combined  acoustic  and  flight  loads.  The  crack  branching, 
however,  is  attributed  to  the  simulated  acoustic  loading  on  comparing  the 
straight  crack  in  test  specimen  CF7  (Figures  57  and  71 ) with  the  cracks  in 
test  specimens  CF5 , CF9-2,  and  CF11  (Figures  58,  60,  72  and  73).  Crack  branch- 
ing also  occurred  in  test  specimens  subjected  only  to  acoustic  loading,  as 
discussed  in  Reference  16. 

The  crack  growth  in  all  the  test  specimens  subjected  to  the  lower  wing  surface 
loading  exceeded  the  calculated  crack  growth  without  retardation,  as  illus- 
trated in  Figure  77  for  test  specimen  CF6  and  in  Figure  78  for  test  specimens 
CF8,  CF10  and  CF12.  Crack  surface  rotation  to  U5  degrees  relative  to  the 
panel  surface  is  characteristic  for  these  test  specimens,  except  for  CF6, 
associated  with  the  relatively  high  flight-by-flight  loads  as  summarized  in 
Table  11.  The  accelerated  crack  growth  and  the  rotation  of  the  crack  face  are 
directly  associated  with  the  high  panel  stresses  (See  Table  11 ) produced  by 
the  high  lower  wing  surface  flight  loads.  The  rotation  of  the  crack  face 
indicates  a change  over  from  a plane  stress/plane  strain  condition  to  a plane 
stress  condition  through  the  thickness  of  the  panel  ahead  of  the  crack  tip 
(Reference  1+0).  Also,  no  crack  branching  was  observed  leading  to  the  conclusion 
that  the  relative  magnitude  of  the  acoustic  load  stress  to  the  tension  flight 
load  stress  is  a significant  parameter  in  crack  branching. 

The  limited  crack  growth  data  from  the  simulated  acoustic  loading  did  indicate 
a considerably  lower  than  expected  crack  growth  rate  as  discussed  in  Sec- 
tion V 3.6.3  and  illustrated  in  Figure  70.  The  reduced  crack  growth  could 
indicate  that  the  expression  for  the  Mode  I stress  intensity  factor  (equa- 
tion 66),  for  a crack  along  a fixed  boundary,  may  not  be  applicable  to  riveted 
panels  where  the  peak  bending  moment  occurs  alternately  at  the  rivet  line  and 
the  frame  edge  in  each  stress  cycle. 
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VII.  CONCLUSIONS 


A.  The  electromagnetic  shaker  arrangement  used  in  this  program  provides  in 
general,  a good  simulation  of  the  acoustic  loading  for  the  fundamental 
mode  of  a stiffened  panel. 

B.  The  test  specimen  design  objectives  with  regard  to  critical  panel 
buckling  stress,  uniform  stress  distribution  across  the  panel,  plane 
stress  state  and  small  variability  between  test  specimens  have  been 
accomplished. 

C.  Good  agreement  between  the  measured  baseline  constant  amplitude  crack 
growth  data  and  comparable  published  data  for  the  same  load  conditions 
has  been  obtained. 

D.  The  axially  applied  flight  loads  produced  a plane  stress  field  in  the 
test  specimen  which  is  parallel  to  the  panel  edge  and  normal  to  the 
rivet  hole  precrack.  The  acoustic  loading  produced  a reverse  bending 
type  stress  distribution  in  the  panel  edge  in  the  vicinity  of  the  pre- 
crack, with  the  peak  surface  stress  path  normal  to  the  panel  edge  and 
the  flight  load  stress  field. 

E.  For  the  structural  model,  crack  orientation  and  combination  of  acoustic 
and  flight  loads  used  in  this  program,  the  measured  crack  growth  data 
did  not  indicate  any  significant  effects  on  crack  growth  by  the  inclu- 
sion of  the  acoustic  loading.  In  other  structures  in  which  the  peak 
acoustic  load  stress  is  acting  in  the  same  direction  as  the  flight  load 
induced  stress  and  the  precrack  is  oriented  normal  to  both  stress  paths, 
significant  effects  on  crack  growth  may  occur. 

F.  Surface  crack  branching  is,  however,  caused  by  acoustic  loading  and  is 
a function  of  the  relative  acoustic  and  flight  tension  load  induced 
stress  levels. 

G.  The  measured  crack  growth  rate  data  on  panels  subjected  to  the  upper 
wing  surface  flight  loads  are  in  good  agreement  with  the  predicted 
crack  growth  rate. 

H.  The  higher  level  lower  wing  surface  loads  produce  crack  face  rotation 
relative  to  the  panel  surface  associated  with  a plane  stress  condition 
ahead  of  the  crack  tip.  The  measured  crack  growth  rate  exceeded  the 
predicted  crack  growth  rate  without  retardation. 

I.  The  limited  acoustic  crack  growth  data,  from  the  acoustic  fatigue  type 
rivet  hole  precrack  parallel  to  the  panel  edge,  indicated  a lower  than 
expected  crack  growth  rate.  Further  evaluation  is  prevented  by  the 
lack  of  random  acoustic  crack  growth  data  for  U.S.  airplane  alloys  and 
for  riveted  airplane  structures,  in  general. 
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